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We  present  here  a  detailed  theoretical  and  experimental  investigation 
of  the  dynamics  of  surface  melting  and  regrowth  and  electrical  activation 
of  Si  and  BF2  ion  implanted  amorphized  silicon  annealed  with  a  Q-switched 
Nd : glass  laser  (X  ■  1.06  ym) .  In  the  theoretical  calculations  we  have  used 
the  technique  of  finite  difference  equations  to  solve  the  one-dimensional 
heat  conduction  equations  for  the  solid  and  liquid  phases.  We  have  taken 
into  account  the  temperature  dependence  of  all  the  parameters  involved  to 
give  us  a  realistic  model.  The  melting  model  is  supported  by  time  re¬ 
solved  reflectivity  measurements  using  a  He-Ne  and  an  Ar  laser.  We  use  SEM 
to  investigate  the  electron  channelling  pattern  to  study  the  quality  of 
crystalline  regrowth.  For  BF2+  implanted  silicon,  the  boron  and  fluorine 
atomic  profiles  are  measured  by  secondary  ion  mass  spectrometry  (SIMS)  on 
as-implanted  samples  as  well  as  laser  annealed  samples.  The  impurity  re¬ 
distribution  is  explained  on  the  basis  of  the  melting  model.  Differential 
resistivity  and  Hall  effect  measurements  in  conjunction  with  successive 
layer  removal  are  used  to  obtain  the  electrical  carrier  distributions  as 
a  function  of  the  laser  fluence.  For  good  quality  crystalline  regrowth 
and  full  electrical  activation,  the  laser  fluence  should  lie  between  the 
melting  threshold  and  the  damage  threshold.  Therefore,  a  reliable  model 
is  useful  in  designing  experiments  of  laser  annealing. 

Pulsed  C02  lasers  (X  »  10.6  ym)  have  the  distinct  advantage  of  having 


considerably  larger  spot-size  than  a  Nd: glass,  YAG  or  Ruby  Laser,  and  are 
therefore  very  promising  for  laser  annealing  of  large  Si  wafers  commonly 
used  in  Industrial  production.  At  10.6  pm,  the  absorption  coefficient  of 
amorphous  silicon  is  much  smaller  than  that  of  crystalline  silicon  having 
large  donor  concentration.  Hence,  for  annealing  with  a  CO^  laser  it  is 
advantageous  to  implant  the  sample  with  a  moderately  high  ion  dose  without 
driving  the  sample  amorphous. 

We  present  a  calculation  of  the  temperature  rise  and  investigate  the 
"thermal  runaway"  phenomenon  during  pulsed  CO^  laser  annealing  of  silicon. 
The  calculations  are  based  on  the  thermal  melting  model,  taking  into 
account  the  temperature  dependence  of  all  pertinent  material  parameters, 
including  the  absorption  coefficient.  In  calculating  the  temperature 
variation  of  free  carrier  absorption  in  n-Si,  we  have  taken  into  account 
acoustic  deformation  potential  scattering,  optical  deformation  potential 
scattering,  and  ionized  impurity  scattering.  The  deformation  potentials 
are  adjusted  to  fit  the  experimentally  observed  values  at  300°K. 
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I.  INTRODUCTION 

The  technique  of  ion  implantation  In  semiconductors  for  device 
fabrication  is  accompanied  by  radiation  damage.  Conventionally,  thermal 
annealing  has  been  employed  to  restore  the  lattice  to  crystalline  order 
and  to  provide  for  electrical  activation  of  the  dopant  atoms.  It  has  been 
shown  recently  [1]  that  an  intense  laser  light  of  short  pulse  duration  has 
several  advantages  over  thermal  annealing  of  ion-implanted  amorphlzed  thin 
layers  in  Si.  We  present  below  a  detailed  account  of  a  thermal  melting 
model  to  explain  the  physical  mechanism  behind  pulsed  laser  annealing. 

The  model  is  experimentally  supported  by  time-dependent  optical  reflec¬ 
tivity  measurements.  Similar  calculations  have  been  done  with  slightly 
different  assumptions  for  ion-implanted  silicon  [2]  and  crystalline 
silicon  [3] .  The  present  calculations  are  compared  with  results  of  time- 

resolved  reflectivity  measurements  for  pulsed  annealing  at  1.06  ym  of  an 

+  *f 

amorphous  layer  produced  by  Implants  of  Si  in  Si  and  BF^  In  Si.  Com¬ 
parison  of  the  calculations  and  the  experimental  model  support  the 
assumption  that  the  Si  surface  melts  during  pulsed  laser  annealing  of 
appropriate  power  density.  The  thermal  melting  model  is  also  supported 
by  atomic  Impurity  profiling  using  the  SIMS  technique,  described  by  Wang 
et  al.  [4]  and  analysis  of  electron-beam  Induced  current  using  SEM,  de¬ 
scribed  by  Jastrzebskl  et  al  [  5]  . 

For  good  quality  crystalline  regrowth,  the  laser  fluence  should 
lie  between  the  melting  threshold  and  the  damage  threshold.  It  is  also 
Important  to  adjust  the  laser  r  'rd  the  mat  rial  parameters  such  that  the 


melt  front  propagates  through  the  amorphous-crystalline  interface  and  wets 
the  substrate.  Otherwise,  polycrystalline  growth  occurs.  Therefore,  a 
predictive  model  is  very  useful  in  choosing  appropriate  experimental  condi¬ 
tions. 

For  a  particular  thickness  of  the  amorphous  layer,  the  incident 
laser  energy  and  the  pulse  width  of  the  laser  beam  are  used  as  parameters 
to  determine  the  duration  of  the  surface  melting  and  the  time-evolution  of 
the  melt  front  into  the  material.  The  numerical  technique  used  in  solving 
the  nonlinear  heat  conduction  equations  is  the  finite  difference  method. 

The  formalism  can  be  applied  to  Ge,  GaAs,  and  other  materials  and  also  to 
various  lasers  such  as  ruby,  CC^j  etc.  The  material  parameters  and  the 
absorption  mechanisms  for  the  different  laser  wavelengths  must  be  taken 
into  consideration,  but  the  theoretical  methods  are  similar. 

We  note  that  in  the  visible  and  near  infrared  region,  the  absorp¬ 
tion  coefficient  of  amorphous  silicon  [6]  is  an  order  of  magniture  or  more 
higher  than  that  of  crystalline  silicon.  We  take  this  factor  into  consid¬ 
eration  in  determining  the  threshold  energy  density  of  melting  of  ion-im¬ 
planted  amorphized  silicon. 

The  phenomenon  of  pulsed  laser  annealing  and  the  experimental 
results  of  time-resolved  reflectivity  measurements  are  explained  on  the 
basis  of  a  melting  model.  There  have  been  speculations,  by  Van  Vechten 
et  al.  [7]  from  theoretical  considerations  and  by  Lo  and  Compaan  [8]  using 
Raman  measurements,  about  the  possible  role  of  a  dense  electron-hole  plasma 
which  might  prevent  the  lattice  temperature  from  rising  significantly 
during  pulsed  laser  annealing  of  silicon.  ,Towever,  we  find  good  agreement 
between  the  rather  simple  melting  model  and  results  of  reflectivity  experi- 


ments,  without  invoking  plasma  effects.  The  present  work  and  time-resolved 
reflectivity  measurements  obtained  by  Nathan  et  al. ,  [9]  suggest  that  the 
rise  in  reflectivity  of  the  Si  surface  and  the  resulting  recrystallization 
is  caused  by  thermal  melting  for  typical  pulsed  laser  annealing  experiments, 
and  is  most  likely  not  a  result  of  the  formation  of  a  dense  electron-hole 
plasma  which  prevents  a  significant  rise  in  temperature. 

1.1  Melting  Model 

We  begin  with  the  ion  implanted  amorphized  silicon  at  room 
temperature  in  the  solid  phase.  The  laser  light  is  then  absorbed  in  the 
near  surface  region  of  the  solid.  The  power  absorbed  by  the  solid  can  be 
expressed  by  a  source  term  which  is  both  spatial  and  time  dependent.  The 
absorbed  light  energy  is  then  converted  into  heat,  which  diffuses  into  the 
bulk  by  thermal  conduction.  The  above  process  can  be  explained  by  consider¬ 
ing  the  heat  conduction  equation  with  an  appropriate  source  term. 

The  temperature  at  the  surface  rises  and  eventually  reaches  the 
melting  point  of  silicon.  The  melt  front  propagates  into  the  material, 
reaches  a  maximum  value  and  then  recedes  toward  the  surface,  leading  to  a 
liquid-phase  epitaxial  regrowth.  When  we  have  both  the  liquid  and  solid 
phases  present,  we  have  to  solve  simultaneously  two  heat  conduction  equations 
for  the  two  phases.  The  position  of  the  melt  front  is  also  calculated  at 
every  time  step.  Finally,  we  have  crystalline  silicon  in  the  solid  phase 
at  room  temperature. 

Since  the  typical  laser  spot  is  about  3  mm  in  diameter  and  the 
pulse  is  uniform,  it  is  appropriate  to  neglect  the  radial  flow  of  heat  and 
consider  only  the  one-dimensional  flow  of  heat  into  the  material. 


The  one-dimensional  heat  flow  is  described  by  the  equation 

C(T)  pit "  Ek<t)  H ]+  s(x»t)  (1-1) 

where  <  is  the  thermal  conductivity,  C  is  the  specific  heat  capacity,  and 
p  is  the  density  of  the  material.  The  source  term  S  is  governed  by  the 
mechanism  of  the  laser-solid  interaction  and  the  time  shape  of  the  laser 
pulse.  The  laser  pulse  is  assumed  to  be  Gaussian,  and  so  the  time-dependent 
part  of  the  source  term  is  given  by 


S(t) 


( 


t-t* 

r 


2 
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(1.2) 


where  T  is  the  full  width  duration  of  the  pulse  at  1/e  intensity  and  t' 
is  the  time  when  the  pulse  reaches  its  maximum  value. 

After  the  front  surface  has  melted,  the  heat  flow  into  the  liquid 
is  described  by  an  equation  similar  to  (1.1)  with  the  appropriate  parameters 
for  the  liquid  phase. 

The  boundary  conditions  employed  are 
(a)  as  x  00 ,  T  300°K 


3T 

(b)  *  0  at  x  »  0.  The  boundary  condition  must  be  appropriately 

modified  to  consider  evaporation  effects  at  the  highest  laser  energies. 

Since  the  thermal  conductivity  and  specific  heat  capacity  are  both  tempera¬ 
ture  dependent,  equation  (1.1)  is  essentially  nonlinear.  We  employ  the 
technique  of  finite-difference  equations  to  solve  equation  (1.1)  for  both 
the  solid  and  liquid  phases,  taking  into  account  appropriate  time-dependent 
boudary  conditions  to  find  the  position  of  the  melt  front  (xffl)  as  a  function 


of  time. 
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At  any  time,  x^  is  obtained  from  solution  of  the  equation 


3T 
^s  3x 


3T 
KL  3x 


+  L 


dx 
_ m 

dt 


(1.3) 


where  k  ,  <T  are  the  thermal  conductivities  of  the  solid  and  liquid  phases 
and  L  is  the  latent  heat  of  melting.  Equation  (1.3)  was  used  in  the 
numerical  computations  to  study  the  time  evolution  of  the  melt-front. 


We  note  that  we  are  considering  a  first-order  phase  transition 
involving  latent  heat  with  the  density  p  as  the  order  parameter.  Hence  the 
density  p  cannot  be  considered  to  be  constant  during  the  phase  change  as 
was  assumed  by  Schultz  and  Collins  [10  ] . 

We  have  considered  the  temperature  variation  of  thermal  conduc¬ 
tivity  [11]  and  specific  heat  capacity  [12]  of  Si  in  the  solid  phase.  The 

variation  of  k  with  temperature  turns  out  to  be  crucial  in  calculating 
s 

the  temperature  rise  during  laser  annealing.  The  thermal  conductivity  of 
liquid  Si  has  the  value  0.322  W/cm-°K  [13]  .  Surko  et  al.  [3]  used  a  value 
of  obtained  from  using  the  Wiedemann-Franz  law.  However,  in  view 
of  the  small  proportion  of  the  electronic  contribution  for  the  near  melting 
solid,  any  estimation  of  the  thermal  conductivity  of  molten  silicon  from 
the  electrical  conductivity  is  not  a  good  approximation.  In  our  work  the 
specific  heat  capacity  of  liquid  Si  was  assumed  [12]  to  have  a  constant 
value  of  0.915  j/g-°K,  since  its  variation  with  temperature  is  not  signifi¬ 
cant.  In  our  calculations  we  have  used  the  values  for  crystalline  silicon 
when  considering  the  temperature  variation  of  the  thermal  conductivity 


and  specific  heat  capacity  of  the  ion- Implanted  amorphlzed  silicon  layer. 

There  are  no  available  data  on  the  values  of  k  and  C  for  amorphous 

s  s 


silicon  from  300°K  to  1685°K.  However,  the  difference  in  the  thermal 

properties  between  crystalline  and  amorphous  silicon  is  less  pronounced  at 

room  temperature  than  at  low  temperatures,  and  becomes  negligible  as  the 

sample  temperature  rises  from  300°K  to  1685°K.  For  liquid  silicon,  the 

values  of  and  are  independent  of  whether  the  sample  melted  from  an 

amorphous  or  a  crystalline  phase  in  the  solid  state.  The  density  in  the 

3  3 

solid  and  liquid  phases  is  2.329  g/cm  and  2.53  g/cm  ,  respectively.  To 
have  the  calculations  and  experiment  agree  reasonably  well,  the  temperature 
variations  of  all  the  parameters  must  be  considered  fully.  Schultz  and 
Collins  [10]  assumed  the  specific  heat  of  silicon  to  be  constant  over  the 
whole  temperature  range.  Baeri  et  al.  [14]  have  considered  the  variation 
of  the  thf rmal  conductivity  of  silicon  in  the  solid  phase,  but  have 
assumed  the  other  parameters  to  be  constant. 

Before  melting,  the  spatial  part  of  the  source  term  is  given  by 

-a  x 

S  (x)  -  E(l-R  )  a  e  S  (1.4) 

s  s  s 

2 

where  E  is  the  fluence  of  the  laser  pulse  in  J/cm  ,  R  is  the  reflectivity, 

s 

and  ag  ig  the  absorption  coefficient.  Ideally,  one  has  to  consider  the 

effect  of  the  crystalline  substrate  when  the  thickness  of  the  implanted 

amorphized  layer  is  less  than  ct  We  note  that  our  present  assumption 

s 

will  give  values  of  calculated  temperatures  somewhat  higher  than  the 
experimental  situation.  At  1.06  um,  the  value  of  the  absorption  coefficient 
for  crystalline  silicon  at  300°K  is  about  10/cm  and  Increases  considerably 
as  the  sample  temperature  rises  during  pulsed  laser  annealing.  This 
temperature  dependent  absorption  coefficient  introduces  an  additional  non¬ 
linearity  and  consequently  further  complications  in  the  calculation.  One 


can  qualitatively  appreciate  the  fact  that,  as  the  thickness  of  the 
amorphous  layer  increases,  the  effect  of  the  crystalline  substrate  de¬ 
creases.  Finally,  when  the  thickness  of  the  amorphous  layer  is  two  to 

three  times  higher  than  a  one  can  possibly  neglect  the  substrate 

s 

effect  without  much  loss  of  accuracy. 

After  melting,  the  source  term  has  an  expression  sililar  to  Eq. 
(1.4),  with  as  the  reflectivity  and  the  absorption  coefficient  of 
liquid  silicon. 


SL(x)  »  E(l-R^)  c^e 


(1.5) 


For  the  case  of  a  Nd: glass  laser,  the  values  of  Rg  and  R^  were 

taken  to  be  0.31  and  0.76,  respectively  [15,16].  This  large  change  in  the 

value  of  reflectivity  during  phase  transitions  was  considered  throughout 

the  calculations  and  was  monitored  experimentally  during  time-resolved 

reflectivity  measurements.  At  1.06  um,  the  values  of  ag  [15]  for  a-Si 

3  6 

and  [16]  for  liquid  silicon  are  given  by  6  x  10  /cm  and  1.4  x  10  / 

cm,  respectively.  We  note  that  the  value  of  is  about  two  orders  of 

magnitude  higher  than  a  .  Therefore,  after  the  front  surface  has  melted, 

s 

we  neglect  the  transmission  of  laser  energy  through  the  liquid  phase. 

This  is  justifiable  as  the  value  of  ^  is  about  70  X  and  all  the  laser 
energy  is  absorbed  by  the  thin  liquid  layer. 

To  estimate  the  effect  of  transmission  of  the  incident  laser 
energy  from  the  liquid  to  the  solid  phase,  one  can  neglect  changes  in  the 
reflection  coefficient  during  melting  and  write  for  the  source  terms  in  the 
solid  (s)  and  liquid  (L) : 


S  (x)  -  E  (1-R  )  a  e 


(1.6) 


V*>  ’  V1-*.)  °Le 


(1.7) 


To  determine  Eg  and  E^  in  the  spirit  of  the  above  approximation, 
we  apply  the  following  boundary  conditions: 

(a)  At  the  junction  of  separation  of  the  two  phases,  the  source 


term  should  be  continuous. 


at  x  -  x  ,  S  (x  )  -  ST (x  ) 
m’  s  m  L'  m 


(1.8) 


(b)  The  absorbed  total  energy  density  E  is  distributed  between 
the  liquid  and  solid  phases  (and  does  not  change  during  melting) , 


E  *  /  "  S.  (x)  dx  +  J  S  (x)dx 

O  L  X  S 

m 


(1-9) 


From  equations  (1.6  -  1.9)  we  obtain: 


-a  (x-x) 
s  m 


Ss(x) 


E(l-R  )  a  e 

s  s _ 

a  OLX 

1  +  —  (e^  m-l) 


(1.10) 


-OL  (x-x  ) 

S.M  .  E(1-V 

a  cl.  x 
1  +  —  (e^*  m-l) 

°L 


(1.11) 


We  point  out  that  using  (1.10)  and  (1.11)  would  prevent  taking 
into  account  the  important  difference  betweeen  Rg  and  R^.  Equations 
(1.10)  and  (1.11)  were  not  Included  in  the  actual  computations.  Since  the 
large  difference  between  Rg  and  R^  is  a  major  concept  in  the  melting  model 
we  have  used  equations  (1.4)  and  (1.5)  in  our  computations.  Bell  [2]  has 


9 


adopted  a  different  approach  in  considering  the  transmission  of  laser  energy 
from  the  liquid  to  solid.  When  the  liquid  layer  exceeds  a  value  of  300 
Bell  assumes  that  all  of  the  incident  light  that  is  not  reflected  is  absorbed 
within  the  first  300  X-thick  meshpoint.  For  a  liquid  layer  of  thickness 
less  than  300  X  the  details  of  transmittance  is  taken  into  account.  How¬ 
ever,  as  pointed  out  by  Bell,  this  detailed  consideration  during  initiation 
of  the  melted  layer  does  not  have  a  large  effect  on  the  results  of  the 
computation.  Since  for  a  laser  fluence  of  2.0J/cm  the  melt-front  reaches 
a  maximum  value  of  more  than  4000  X,  the  details  of  the  transmission  in  the 
first  300  X  liquid  layer  should  be  rather  Insignificant.  As  discussed 
above,  since  the  value  of  (1.4  x  10^/cm  )  is  much  larger  than  the  value 

3 

of  ag  (6  x  10  /cm  ),  we  are  justified  in  neglecting  the  transmission  of 
laser  energy  from  liquid  to  solid  after  the  front  surface  has  melted. 

To  solve  equation  (1.1)  numerically  we  divide  the  sample  Into 
equal  sections  of  length  Ax.  The  time  step  is  denoted  by  At.  The  intervals 
Ax  and  At  are  chosen  to  satisfy  the  von  Neumann  stability  criterion  given 
by 

At  <  —  (Ax)2 

Specifically,  we  choose  Ax  ■  250  X  and  At  »  2  psec  for  our  calcu¬ 
lation.  The  thickness  of  a  sample  was  chosen  to  be  50  ym  for  numerical 
computations.  If  n  and  m  denote  a  particular  set  of  spatial  and  time 
coordinates,  then  a  typical  expression  for  the  temperature  as  a  solution 
to  equation  (1.1)  is  given  by 


t- 


m 
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m 


1“  -  t“-1  + 

n  n 


±  ( 
Ax  \ 


t”:,1  - 1”-1  r-1  -  T" 

n+1  n  n 

sl  Ax  s2  Ax 


^-1  V 

■=0 


(1.12) 


+  S  (x  ,t  ) 
s  n*  m' 


At 


p  C  (t“_1) 
s  s  n 


Similar  equations  govern  the  temperature  rise  in  the  liquid  phase. 
1. 2  Ion  Implantation 

Ion  implantation  is  now  a  standard  technique  of  doping  semicon¬ 
ductors  and  has  been  extensively  reviewed  in  the  literature  [17-21]. 
Basically,  ion  implantation  consists  of  energetic  ions  entering  the  sub¬ 
strate  and  losing  energy  through  collisions  with  the  lattice  until  they 
come  to  rest  in  the  substrate.  The  distribution  of  ions  in  the  substrate 
is  a  function  of  the  incident  ion Species  and  its  energy,  along  with  prop¬ 
erties  of  the  substrate  material.  The  two  major  forms  of  energy  loss  for 
the  incoming  ion  are  interactions  with  the  electrons  of  the  substrate  atoms 
and  collisions  with  the  nuclei.  While  only  approximately  true,  these  two 
loss  mechanisms  can  be  assumed  to  be  independent  of  one  another.  Under 
this  assumption  we  can  write 


-  N  [sn<E>  +  Se(E)j  (1.13) 

where  E  is  the  energy  of  the  ion  at  a  point  x  along  its  path,  Sq(E)  is  the 
nuclear  stopping  power.  S^CE)  is  the  electronic  stopping  power,  and  N  is 
the  average  number  of  target  atoms  per  unit  volume. 


w 


To  a  fairly  good  approximation,  the  nuclear  stopping  power  can 


be  assumed  to  be  elastic  scattering  processes,  independent  of  one  another. 

In  the  context  of  these  assumptions,  the  collisions  can  be  treated  with 
classical  mechanics.  When  a  potential  which  takes  into  account  the  screening 
of  electrons  is  used,  fairly  accurate  values  of  the  nuclear  stopping  power 
can  be  obtained. 

The  electronic  stopping  power  is  more  difficult  to  compute 
theoretically.  There  are  basically  two  different  models  that  have  been 
used  to  calculate  this  quantity.  The  first  method  approximates  the  elec¬ 
trons  in  the  target  by  a  free  electron  gas.  In  this  approximation  the  ion 
loses  energy  by  collisions  with  the  electrons  and  excitation  of  plasma 
resonances.  For  ion  velocities  less  than  that  of  an  electron  having  an 
energy  below  the  Fermi  energy  of  the  gas,  Lendhard  and  Winther  [22]  have 
shown  that  Se(E)  can  be  given  by 

S  (E)  *  CV  «  kE5  (1.14) 

e 


where  the  parameter  k  is  a  function  of  the  ion  and  substrate  involved. 

The  other  model,  proposed  by  Firsov  [23  ]  ,  is  based  on  the  elec¬ 
tronic  interactions  in  a  two-body  collision.  The  process  may  be  visualized 
by  the  two  atoms  forming  a  "quasi-molecule"  as  they  approach  each  other. 

The  atoms  will  exchange  outer  electrons  as  they  attempt  to  form  a  molecule. 
These  electrons  will  take  on  the  net  momentum  of  the  "quasi-molecule". 

When  the  atoms  separate  again,  the  incoming  ion  will  have  lost  energy  to 
the  stationary  one.  While  the  Firsov  model  assumes  participation  of  all 
the  electrons,  it  is  more  appropriate  to  assume  some  screening  from  the 


inner  electrons.  Electronic  stopping  power  is  more  important  for  high 
energies  and  light  ions,  while  nuclear  stopping  is  the  major  force  for  low 
energies  and  heavy  ions. 


Once  Sq  and  Sg  are  known,  they  can  be  substituted  into  equation 
1.13,  and  the  total  distance  or  range  (R)  the  ion  travels  on  its  winding 
path  in  the  substrate  can  be  found.  Experimentally,  we  can  only  observe 
the  distance  travelled  perpendicular  to  the  surface,  the  projected  range 
(Rp) .  Also,  since  the  ions  suffer  collisions  in  a  random  fashion,  the 
ions  will  come  to  rest  in  some  type  of  distribution  about  the  projected 
range.  The  spread  of  the  distribution  is  defined  by  the  standard  deviation 
(AR) ,  and  the  spread  perpendicular  to  the  surface  is  described  by  the  pro¬ 
jected  standard  deviation  (AR^) . 

Differential  equations  for  computing  Rp  and  AR^  have  been  derived 
by  Lindhard,  Scharff ,  and  Schiott  [24]  .  This  theory,  henceforth  referred 
to  as  LSS,  describes  the  basic  features  of  the  implanted  profile.  The 
Implanted  profile  can  usually  be  approximated  by  a  gaussian  function  given 
by 


N(x)  ■  Np  exp 


■<»  -  V' 

2AR 


(1.15) 


_3 

where  N(x)  is  the  concentration  of  atoms  (cm  )  at  a  point  x  and  Np  is  the 
peak  concentration.  This  peak  concentration  is  expressed  by 

0.4 


N  .  N 
p  AR  D 
P 


(1.16) 


where  is  the  fluence  or  dose  of  the  implant  expressed  in  atoms/cm  . 
Values  of  Rp,  ARp,  and  related  quantities  have  been  tabulated  for  various 
ions  and  substrates  [25].  Distortion  of  the  profile  when  higher  order 
moments  are  considered  has  also  been  described  [25] . 
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While  the  LSS  theory  gives  fairly  accurate  impurity  profiles  for 
unannealed  samples,  the  distribution  can  change  drastically  after  the 
annealing  process  necessary  to  restructure  the  crystal  and  activate  the 
Impurities.  Considerable  work  has  been  done  to  study  how  various  impuri¬ 
ties  diffuse  during  annealing  and  how  well  the  impurities  are  electrically 
and  optically  activated.  Both  of  these  pieces  of  information  are  essential 
in  fabricating  devices  with  the  desired  doping  profiles. 

The  implantation  process  consists  basically  of  forming  a  plasma 
containing  the  desired  ions,  accelerating  the  ions  to  a  predetermined 
energy,  selecting  only  the  appropriate  ion  species  out  of  the  beam  by  mass 
separation,  and  directing  the  ions  uniformly  onto  the  substrate  by  scanning 
the  beam  in  some  manner.  The  charge  hitting  the  sample  is  countered  so 
that  the  desired  dose  of  impurities  can  be  obtained.  A.  schematic  diagram 
of  the  300  keV  Accelerators,  Inc.  300-MP  ion  lmplanter  used  in  this  work 
is  shown  in  Figure  1.1. 


Figure  1.1.  Schematic  diagram  of  the  ion  implantation  system. 


2.  LASER  ANNEALING  OF  Si+  IMPLANTED  AMORPHIZED  SILICON 

In  this  section  we  study  theoretically  and  experimentally  the 
dynamics  and  quality  of  crystalline  regrowth  during  annealing  of  multiple 
Si+  ion  implanted  amorphized  silicon  with  a  Q-switched  Nd:glass  laser. 

Crystalline  silicon  can  be  transformed  into  amorphous  Si  by 
appropriate  high-dose  Implantation  [26]  .  There  have  been  two  models  pro¬ 
posed  to  calculate  the  critical  doses  for  this  transformation:  a  critical- 
energy-density  (CED)  model  [27,28  ],  and  a  damage-overlap  (DO)  model 
[29,30].  Since  amorphous  layers  play  an  Important  role  in  annealing  be¬ 
havior,  the  prediction  of  formation  of  amorphous  Si  is  very  Important,  as 
is  it's  detection  after  Implantation  [31,32]. 

We  first  present  theoretical  results  based  on  the  thermal  melting 
model  for  ion  implanted  amorphized  silicon. 

We  describe  the  experimental  set  up  of  the  Nd: glass  laser  system 
for  annealing  experiments. 

The  experimental  set  up  for  time  resolved  reflectivity  measurements 
and  the  experimental  results  are  then  presented. 

To  study  the  quality  of  the  crystalline  regrowth  we  present  a 
electron  channeling  pattern  using  a  Scanning  Electron  Microscope. 

2.1  Theoretical  Results 

A  typical  temperature  vs.  time  profile  for  self-ion  implanted 

amorphized  silicon  irradiated  with  a  Q-switched  Nd: glass  laser  of  50  nsec 

2 

FWHM  and  laser  fluence  1.5  J/cm  is  shown  in  Fig.  2.1.  The  shape  of  the 

laser  pulse  is  shown  in  the  bottom  left-hand  corner.  T  denotes  the 

in 

melting  point  of  silicon,  which  is  1685°K.  With  the  laser  pulse,  we  find 
that  the  temperature  at  the  surface  rises  rapidly  with  a  parabolic  trend 
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Figure  2.1.  Results  of  the  calculation  showing  the  temperature  vs.  time 

profiles  at  the  surface,  4000  A,  and  8000  A  layers  for  ion- 

implanted  amorphized  silicon  irradiated  with  a  Q-switched 

2 

Nd: glass  laser  of  fluence  1.5  J/cm  and  50  nsecs  FWHM. 


until  the  surface  melts.  This  is  consistent  with  the  simple  phase  change 

model  as  discussed  by  Carslaw  and  Jaeger  [33]  .  The  surface  melts  after 

the  laser  pulse  reaches  its  peak.  After  melting,  the  temperature  of  the 

liquid  rises  rather  sharply.  The  temperature  profiles  at  depths  4000  £ 

2 

are  also  plotted.  We  note  that  with  the  laser  fluence  of  1.5  J/cm  ' ,  the 

4000  2.  and  8000  8.  layers  never  melt.  Therefore,  to  anneal  a  4000  % 

2 

amorphous  layer  of  Si,  a  laser  fluence  of  1.5  J/cm  is  not  sufficient. 

This  is  confirmed  by  experiment.  We  also  note  that  for  deeper  levels,  the 

maximum  temperature  occurs  at  later  times,  which  is  expected  from  the 

mechanism  of  thermal  diffusion  of  heat  into  the  material. 

In  Fig.  2.2  we  show  a  similar  plot  of  temperature  vs.  time  for 
2 

a  laser  fluence  of  2.0  J/cm  .  In  this  case,  the  surface  melts  at  about 
the  peak  of  the  laser  pulse.  Therefore,  as  the  laser  fluence  is  increased, 
the  initiation  of  melting  moves  toward  the  leading  edge  of  the  pulse. 

Also,  the  temperature  of  the  liquid  rises  sharply  to  a  higher  value  and  the 
surface  remains  molten  for  a  longer  period  of  time.  We  also  note  that  in 
this  case  the  4000  %.  layer  melts,  so  that  a  liquid-phase  epitaxial  growth 
of  4000  %.  thickness  of  amorphous  silicon  can  be  expected.  This  example 
illustrates  the  relevance  of  the  model  for  annealing  ion- Implanted 
silicon. 

Another  important  result  of  the  calculation  is  the  time  evolution 
of  the  melt  front  into  the  material  for  a  particular  laser  fluence.  Fig. 
2.3  shows  the  position  of  the  melt  front  x^  vs.  time  for  three  different 
laser  fluences.  We  note  that  as  the  laser  fluence  is  increased,  the  ini¬ 
tiation  of  melting  occurs  earlier,  as  expected.  Also,  the  melt  front 
propagates  deeper  into  the  material  with  increase  of  laser  fluence.  In 
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Figure  2.2. 


Results  of  the  calculation  showing  the  temperature  vs.  time 

2 

profiles  at  different  depths  for  a  laser  fluence  of  2.0  J/cm 
(50  nsecs  FWHM) . 
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particular,  with  a  laser  fluence  of  2.0  J/cm  ,  the  melt  front  propagates 

to  a  maximum  depth  of  about  4000  X  into  the  material.  The  time  at  which 

the  melt  front  reaches  its  maximum  value  shifts  toward  the  decreasing  edge 

of  the  laser  pulse  as  the  laser  fluence  is  increased. 

From  Fig.  2.3  we  can  estimate  the  velocity  of  propagation  (v  ) 

m 

of  the  melt  front  into  the  material  and  also  the  velocity  of  crystalliza¬ 
tion  (v  )  as  the  melt  front  recedes  toward  the  surface.  We  find  that  for 
c 

2 

values  of  laser  fluence  of  1.5,  1.75,  and  2.0  J/cm  ,  the  value  of  v  is 

m 

830,  1250,  1600  cm/ sec  respectively.  We  note  that  with  the  increase  of 
the  laser  fluence,  the  velocity  of  the  melt  front  also  increases.  The 

2 

velocity  of  crystallization  (v^)  for  the  laser  fluences  of  1.5,  2.0  J/cm 

are  found  to  225,  233,  and  240  cm/sec.  respectively.  We  observe  experi- 

2 

mentally  that  for  a  laser  fluence  of  3.0  J/cm  or  more  there  is  visible 
surface  damage. 

The  time  during  which  the  surface  remains  molten  for  a  particular 

value  of  the  laser  fluence  can  be  obtained  from  curves  such  as  Fig.  2.3. 

2 

For  laser  fluences  of  1.5,  1.75,  and  2.0  J/cm  ,  the  duration  of  surface 

melting  (t  )  is  145,  190,  and  235  nsec,  respectively.  These  values  of  t 
m  m 

will  be  compared  below  with  the  experimental  values  obtained  from  the  time- 
resolved  reflectivity  measurements.  We  note  that  the  calculated  velocity 
of  crystallization  depends  sensitively  on  the  assumed  values  of  and  o^. 
The  value  of  0.322  W/cm-°K  used  here  [13]  for  is  less  than  the  value 
of  0.50  W/cm-°K  obtained  by  Surko  et  al.  [3]  from  electrical  resistivity 
data  using  the  Wiedemann -Franz  law.  We  note  that  if  the  assumed  value  of 
k ^  is  greater  than  the  actual  value,  then  the  heat  conduction  into  the 
material  will  be  overestimated.  The  calculated  surface  temperature  will 
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be  too  small  and  the  maximum  value  of  the  melt  depth  will  be  too  large. 

The  increase  in  maximum  value  of  x  will  increase  the  calculated  recrystal- 

m 

lization  velocity.  This  is  demonstrated  by  the  fact  that  calculated  re- 
crystallization  velocity  increases  with  the  maximum  value  of  melt  depth. 
Conversely,  if  the  value  of  is  underestimated,  the  calculated  heat  con¬ 
duction  is  too  small.  As  the  silicon  melts,  the  absorption  coefficient 
increases  rapidly,  since  the  liquid  silicon  behaves  like  a  metal.  Hence, 
assuming  a  constant  average  value  for  the  absorption  coefficient  [14] 
would  lead  to  rather  unjustifiable  values  for  the  velocity  of  crystallization. 
Our  calculations  are  in  principle  similar  to  those  obtained  by 

<3 

Bell  [2] .  Bell,  however,  assumed  constant  values  of  density  p  ■  2.33  g/cmJ 

and  specific  heat  capacity  C  ■  0.95  J/g-°  K  throughout  the  calculations. 

3  3 

We  have  assumed  values  of  density  to  be  2.329  g/cm  and  2.53  g/cm  in  the 
solid  and  liquid  phases  respectively.  Since  we  are  considering  a  first 
order  phase  transition  with  the  density  p  as  the  order  parameter,  the  dif¬ 
ference  in  the  density  between  the  solid  and  liquid  phases  should  be  taken 
into  account.  We  have  also  considered  the  temperature  dependence  of  the 
specific  heat  as  it  varies  from  a  value  of  0.717  j/g-°  K  at  300°  K  to  1.015 
J/g-°  K  below  the  melting  point.  The  specific  heat  capacity  of  liquid 
silicon  is  assumed  to  have  a  constant  value  of  0.915  J/g-°  K.  Other  than 
these  two  parameters  p  and  C,  the  temperature  and  material  dependence  of 
all  the  parameters  have  been  fully  considered  in  our  calculation  and  in 
those  obtained  by  Bell. 

2.2  Sample  Preparation 

To  prepare  self-ion  implanted  amorphized  silicon  we  start  with 
single-crystal  2fl-cm  boron -doped  on  4fi-cm  phosphorus-doped  silicon  samples. 


We  choose  the  energies  and  ion  doses  to  produce  a  continuous  amorphous 
layer  of  ~  4000  X  by  using  a  computer  program  developed  in  our  group  for 
profile  tailoring  [34]  .  The  energies  and  doses  chosen  for  this  work  and 
the  calculated  Si  implant  profile  are  shown  in  Fig.  2.4.  This  series  of 
Si+  implants  was  performed  at  ~  163°  K  to  ensure  amorphization.  For 
these  low  temperature  implants  the  samples  were  attached  to  a  heat  sink 
cooled  by  liquid  nitrogen.  To  reduce  channeling  effects,  all  the  im¬ 
plantations  were  performed  a  few  degrees  ('s'7°)  away  from  the  surface 
normal. 

2.3  Time-Resolved  Reflectivity  Measurements 


In  studying  the  results  of  pulsed  laser  annealing  of  amorphous 
layers,  several  standard  techniques  such  as  Rutherford  backscattering  or 
transmission  electron  microscopy  could  be  used  to  examine  the  quality  of 
the  crystalline  regrowth  mechanism  itself.  Since  the  melting  model  is 
intended  to  describe  the  dynamics  of  pulse  laser  annealing,  it  is  useful 
to  study  the  melting  and  regrowth  as  they  occur.  As  the  silicon  melts 
during  the  laser  pulse,  the  reflectivity  of  the  surface  should  provide  a 
good  monitor  of  the  phase  changes  [15,161.  Therefore,  time-resolved  re¬ 
flectivity  measurements  were  performed  for  comparison  with  the  melting 
model. 

Using  a  photomultiplier  tube  (EMI  9785  B)  of  risetime  1.8  nsec, 
we  have  experimentally  investigated  the  time  dependent  optical  reflectivity 
during  the  surface  melting  process  (Fig.  2.5).  Transitions  from  the 
amorphous  to  crystalline  states  were  achieved  in  Si  with  a  Nd .’glass  pulse 
of  50  nsec  FWHM.  The  combination  of  high  gain,  low  noise,  and  wide  band¬ 
width  makes  the  photomultiplier  particularly  useful  for  studying  fast 


Concentration  (cm"3) 


Figure  2.5.  Experimental  setup  for  time-dependent  optical  reflectivity 


measurements  using  a  He-Ne  laser  during  irradiation  of  the 
self-ion  implanted  amorphized  silicon  with  a  Q-switched 
Nd:glass  laser. 


transients.  Johnson  noise  is  expected  to  be  of  little  consequence  in 
typical  photomultiplier-amplifier  systems. 


We  used  a  5  mW  He-Ne  laser  operating  at  633  nm  and  incident  at 
57.5°  with  the  sample  to  study  the  change  in  reflectance  during  the  annealing 
transient.  The  reflected  light  from  the  surface  was  passed  through  an 
infrared  filter  (3.0.D  at  1.06  ym) ,  a  red  glass  (Schott  G-6300)  to  filter 
out  the  flash  lamps,  and  an  interference  filter  (FWHM  7  nm)  before  reaching 
the  photomultiplier  tube.  Transient  signals  displayed  on  the  oscilloscope 
(Tektronix  7834)  for  a  single  pulse  were  recorded  by  a  Polaroid  camera. 

2.4  Experimental  Results  and  Comparisons  With  the  Melting  Model 
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A  typical  reflectivity  change  at  the  surface  during  a  2.0  J/cm 
Ndrglass  laser  pulse  is  shown  in  the  top  frame  of  Fig.  2.6.  The  laser 
pulse  is  displayed  in  the  bottom  frame.  We  see  that  the  reflectivity 
increases  abruptly  from  its  value  for  amorphous  silicon  at  room  temperature 
(~  40%)  as  the  surface  becomes  hotter.  As  the  laser  pulse  nears  its  peak, 
the  front  surface  melts  and  the  reflectivity  rises  to  a  value  of  65%,  which 
is  consistent  with  the  value  of  reflectivity  for  liquid  silicon  at  633  nm. 
[16] .  The  reflectivity  stays  at  this  value  for  a  time  given  by  the  flat 
portion  of  the  reflectivity  curve.  Then  the  reflectivity  decreases  and 
reaches  a  value  (~  35%)  corresponding  to  that  of  crystalline  silicon  at 
room  temperature.  For  laser  fluences  used  to  observe  melting  of  the  Si 
surface  we  observe  a  fast  rise  followed  by  a  temporary  drop  in  reflectivity. 

This  drop  in  reflectivity  suggests  that  the  reflectivity  of  liquid 
silicon  at  a  higher  temperature  than  the  melting  point  (1685°  K)  is  less 
than  that  of  liquid  silicon  at  the  melting  point.  It  may  be  speculated 
that  this  arises  due  to  change  in  optical  constants  of  liquid  silicon  much 


Figure  2.6.  (a)  Reflectivity  profile  of  the  Si  surface  at  633  nm  as  it 

passes  through  the  phase  transitions  for  a  laser  fluence 
2 

of  2.0  J/cm  . 


(b)  Laser  pulse  (50  nsecs  FWHM) . 


above  the  melting  point.  As  our  calculated  temperature  profile  in  Fig.  2.3 

shows,  there  is  an  initial  sharp  rise  in  the  surface  temperature  in  the 

liquid  phase  after  the  front  surface  has  reached  the  melting  point.  For 

2 

the  laser  fluence  of  2.0J/cm  ,  the  calculated  surface  temperature  does  not 

rise  to  the  boiling  point  (2953°  K) ,  and  there  is  no  experimental  evidence 

2 

of  vaporization  damage  at  2.0J/cm  .  However,  for  laser  fluences  of 
2 

3. OJ/cm  and  above  there  is  visible  surface  damage,  suggesting  that  the 
liquid  silicon  temperature  might  rise  above  the  boiling  point,  leading  to 
vaporization. 

For  various  laser  fluences,  we  measure  the  duration  of  the  flat 
portion  (t  )  of  the  reflectivity  curve,  which  gives  us  the  time  that  the 
surface  remains  molten.  We  can  then  compare  these  values  with  those 
obtained  from  the  plot  of  the  time  evolution  of  the  melt  front  as  given 
by  Fig.  2.3.  The  theoretical  and  experimental  results  are  shown  in  Fig. 

2.7.  We  find  good  agreement  between  the  theoretical  predictions  of  the 
melting  model  and  the  data.  We  also  note  that  the  value  of  the  reflectivity 
of  liquid  silicon  does  not  change  with  the  increase  in  pulse  energy.  How¬ 
ever,  for  laser  fluences  less  than  the  melting  threshold,  the  reflectivity 
does  rise  to  a  value  corresponding  to  that  of  hot  silicon  in  the  solid 
phase. 

The  onset  of  the  reflectivity  change  was  observed  to  move  toward 
the  leading  edge  of  the  pulse  as  the  laser  fluence  was  increased.  This 
is  also  supported  by  the  melting  model  as  shown  by  Fig.  2.3.  After  a 
spot  had  been  recrystallized,  a  second  laser  pulse  of  the  same  fluence  on 
the  same  spot  did  not  produce  any  change  in  reflectivity.  This  is  due  to 
the  fact  that  the  melting  threshold  of  crystalline  silicon  is  considerably 
higher  than  that  of  amorphous  silicon.  At  the  Ndtglass  laser  wavelength 


duration  vs.  laser  fluence  during  laser  irradiation  of  self-ion 
implanted  amorphized  silicon  with  a  Q-switched  Nd: glass  laser 
(50  nsecs  FWHM) . 


the  absorption  coefficient  of  amorphous  silicon  is  about  a  factor  of  100 

more  than  that  of  crystalline  silicon  [6] .  Therefore,  the  amorphous 

sample  absorbs  laser  energy  much  more  efficiently  prior  to  melting  than 

does  the  crystalline  silicon.  This  may  explain  observed  variations  of  t 

m 

with  Implantation  dose  [35]  .  For  laser  fluences  greater  than  the  damage 

2 

threshold  (*“  3.0  J/cm  ) ,  we  find  an  irreversible  decrease  in  reflectivity 
due  to  surface  damage. 

Also  from  Fig.  2.6  we  see  that  after  the  molten  period  (x  )  the 

m 

reflectivity  of  the  surface  decreases  as  the  surface  cools  and  eventually 
reaches  the  value  of  that  of  crystalline  Si  at  room  temperature.  We  note 
that  initially  we  have  amorphous  silicon  of  reflectivity  40%  at  room 
temperature  and  finally  we  have  crystalline  silicon  of  reflectivity  35% 
at  room  temperature.  The  observed  change  in  surface  reflectivity  before 
and  after  the  experiment  is  consistent  with  the  phase  transition  from 
amorphous  to  crystalline  material. 

2.5  Electron  Channeling  Pattern  Using  a  Scanning  Electron  Microscope 

When  a  parallel  electron  beam  rocks  about  a  spot  on  the  specimen 
surface,  the  number  of  backscattered  electrons  varies  according  to  the 
correlation  between  the  incident  beam  and  the  crystal  orientation  of  the 
specimen.  Let  us  assume  that  a  particular  set  of  crystal  planes  satisfies 
the  Bragg  diffraction  condition, 

2d  sin  9  -  nX  (2.1) 

where  d  is  the  interplanar  spacing,  9  is  the  angle  between  the  beam  and  the 
planes,  and  X  the  electron  wavelength. 

In  this  case,  the  penetration  depth  of  the  incident  beam  into 


the  specimen  varies,  and  thus  the  intensity  of  the  backscattered  electrons 


varies  with  respect  to  the  beam  incidence  angle.  Accordingly,  if  the 
incident  beam  is  angle-scanned  in  two  directions,  various  pairs  of  parallel 
lines  referred  to  as  the  Electron  Channeling  Pattern  appear  on  the  cathode 
ray  tube.  These  lines  correspond  to  the  lattice  planes  having  various 
spacings  in  the  crystal. 

To  investigate  the  material  properties  of  the  amorphized  layer 

before  and  after  the  laser  pulse,  we  employed  the  technique  of  electron 

channeling  using  a  scanning  electron  microscope.  For  the  as-implanted 

material  there  was  no  electron  channeling  pattern,  confirming  the  presence 
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of  an  amorphous  layer.  The  sample  which  was  irradiated  with  a  2.5  J/cm 
laser  fluence  clearly  showed  an  electron  channeling  pattern  due  to  the 
crystalline  regrowth  of  the  amorphous  layer.  The  regrown  layer  has  a 
<100>  orientation,  as  expected  for  recrystallization  on  the  <100>  substrate 
A  typical  picture  of  the  channeling  pattern  is  shown  in  Fig.  2.8. 

The  channeling  pattern  shown  in  Fig.  2.8  is  as  sharp  as  that 
obtained  from  an  unimplanted  sample  from  the  same  wafer.  This  supports  the 
fact  that  a  fluence  of  2.5  J/cm  is  sufficient  for  the  melt -front  to  move 
through  the  original  amorphous-crystalline  interface,  as  predicted  from 


the  calculations. 


Figure  2.8.  Electron  channeling  pattern  of  the  Si  implanted 

amorphized  Si  <100>  after  irradiation  with  a  laser 
2 

pulse  of  2.5  J/cm  . 
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3.  LASER  ANNEALING  OF  BF*  IMPLANTED  AMDRPHIZED  SILICON 

Pulsed  laser  annealing  has  been  suggested  as  a  method  to  process 
semiconductor  devices  fabricated  by  ion  implantation.  In  the  fabrication 
of  electronic  devices  in  Si,  the  use  of  high-dose  implantation  is  often 
necessary.  Such  Implantation  produces  an  amorphous  layer  of  a  few  thousand 
angstroms  on  the  silicon  surface.  The  recrystallization  of  such  amorphous 
layers  during  pulsed  laser  annealing  and  the  accompanying  electrical  acti¬ 
vation  are  not  only  of  interest  fundamentally,  but  also  are  important  to 
device  fabrication.  The  physical  mechanism  behind  the  laser-semiconductor 
interaction  leading  to  crystalline  growth  and  electrical  activation  of 
dopant  atoms  is  based  on  a  thermal  melting  model  [36,  2]. 

The  implantation  of  BF2  molecular  ions  into  Si  is  a  useful  method 

for  acceptor  doping,  in  which  B  provides  the  dopant  and  F  accomplishes 

the  amorphlzation  [  31].  The  range  statistics  for  molecular  icn  (BF2) 

implants  can  be  predicted  by  assuming  that  the  BF2  molecule  splits  into 

its  components  with  energies  apportioned  by  mass  ratio.  For  example,  a  150 

19  11 

keV  BF2  molecule  results  in  two  F  with  58  keV  and  one  B  with  34  keV  in 

15  2 

the  Si  sample.  BF2  implanted  Si,  with  fluence  >  1  x  10  /cm  ,  has  high 
electrical  activation  after  low-temperature  thermal  annealing,  and  p+-n 
junctions  formed  by  BFj  implants  have  low  reverse  leakage  current.  For  the 
fabrication  of  devices  requiring  shallow  p  regions,  BF^  is  easier  to 
implant  than  B+,  since  the  implantation  energy  for  BF^  is  much  larger  than 
that  for  B+  to  form  identical  boron  range  distributions,  and  more  stable 
ion  beam  currents  can  therefore  be  obtained  for  BF^  implants.  Another 
advantage  of  BF^  implantation  is  the  amorphization  provided  by  the  heavier 
fluorine  atom.  Conventional  thermal  annealing  at  -  550°C  results  in 
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regrowth  of  the  amorphized  layer,  but  does  not  electrically  activate  the 
tail  of  the  B  distribution,  located  below  the  original  amorphous- 
crystalline  interface  [31].  Only  in  a  thicker  amorphous  layer,  formed  for 
example  by  multiple  Si+  predamage  implants,  can  the  entire  B  profile  be 
activated  [  31  ]. 

We  have  used  the  technique  of  pulsed  laser  annealing  to  electrically 
activate  the  B  distribution  by  liquid  phase  epitaxial  regrowth,  thereby 
avoiding  the  complication  of  implanting  B  into  a  thicker  amorphous  layer 
produced  by  Si+  implants.  The  laser  fluence  corresponding  to  the  melting 
threshold  is  determined  experimentally  by  time-resolved  reflectivity  measure¬ 
ments  [  36,  37].  The  boron  atomic  distributions  are  determined  by  the  SIMS 
technique.  Differential  resistivity  and  Hall  effect  measurements  in  con¬ 
junction  with  successive  layer  removal  are  used  to  obtain  the  electrical 
carrier  distributions  as  a  function  of  the  laser  fluence.  We  find  that  to 
obtain  full  electrical  activation,  the  depth  of  melting  must  exceed  the 
extended  tail  of  the  defect  distribution,  which  is  much  deeper  than  the 
calculated  value.  The  presence  of  such  an  extended  defect  distribution  has 
been  attributed  to  enhanced  diffusion  and  partial  channelling  [  38]. 

3.1  Theoretical  Results 

The  details  of  the  melting  model  have  been  discussed  in  Chapter  1. 

The  laser  light  is  absorbed  in  the  near  surface  region  of  the  solid.  The 
power  absorbed  by  the  solid  can  be  expressed  by  a  source  term  which  is  both 
spatial  and  time  dependent.  The  absorbed  light  energy  is  then  converted 
into  heat,  which  diffuses  into  the  bulk  by  thermal  conduction.  This  process 
can  be  explained  by  considering  the  one-dimensional  heat  conduction  equation 


with  an  appropriate  source  term  given  by 


C(T)pf£  “  <K(T)  +  S(Xjt)  O.l) 

where  <  is  the  thermal  conductivity,  C  is  the  specific  heat  capacity,  and 
p  is  the  density  of  the  material.  The  source  term  S  is  governed  by  the 
mechanism  of  the  laser-solid  interaction  and  the  time  shape  of  the  laser 
pulse.  The  laser  pulse  is  assumed  to  be  Gaussian  with  FWHM  ■  27.5  ns. 

The  source  term  is  given  by 

S(x,t)  -  -r-  exp  E(l-R)  cte"ax  (3.2) 

Tr  r  1 

where  T  is  the  full  width  duration  of  the  pulse  at  1/e  intensity,  t'  is 

the  time  when  the  pulse  reaches  its  maximum  value,  E  is  the  laser  fluence 
2 

in  J/cm  ,  R  is  the  reflection  coefficient,  and  o  is  the  absorption  coef¬ 
ficient.  After  the  front  surface  has  melted,  the  heat  flow  into  the 
liquid  is  described  by  an  equation  similar  to  (3.1),  with  parameters  approp¬ 
riate  for  the  liquid  phase. 

We  have  used  the  technique  of  finite  difference  equations  to  solve 

the  nonlinear  heat  conduction  equations  for  the  solid  and  liquid  phases, 

taking  into  account  appropriate  time-dependent  boundary  conditions  for 

finding  the  position  of  the  melt-front  as  a  function  of  time.  We  have 

considered  the  temperature  variation  of  thermal  conductivity  and  specific 

heat  capacity  of  Si  in  the  solid  and  liquid  phases.  The  density  of  Si  in 

3  3 

the  solid  and  liquid  phases  is  2.33  g/cm  and  2.53  g/cm  ,  respectively. 

For  the  case  of  a  Nd:glass  laser  (X  *  1.06  ym) ,  the  values  of  the  reflection 
coefficient  in  the  solid  and  liquid  phases  are  0.31  and  0.76,  respectively. 
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This  large  change  in  the  value  of  reflectivity  during  phase  transitions 
was  considered  throughout  the  calculations  and  was  monitored  experimentally 
during  time-resolved  reflectivity  measurements. 

At  1.06  ym,  the  value  of  the  absorption  coefficient  of  amorphous  Si 
3 

is  6  x  10  /cm.  Since  the  implanted  amorphous  layer  is  -  1325  A  thick, 
we  have  considered  the  substrate  effects,  especially  the  marked  difference 
between  the  absorption  coefficients  of  amorphous  and  crystalline  Si  at  1.06 
ym.  The  absorption  coefficient  of  the  crystalline  substrate  depends  markedly 
on  temperature  arising  from  the  interband  absorption  and  free-carrier 
absorption  [39].  The  a  for  crystalline  Si  varies  from  14.5/cm  at  300°K 
to  4.7  x  10  /cm  at  the  melting  point  (1685°K).  This  variation  of  a  can 

be  fitted  to  an  empirical  relation  a  (c-Si)  *  1.35  x  10  ^  T^’^/cm  [40]. 

For  a  realistic  model  dealing  with  a  thin  amorphous  layer  of  thickness  much 
smaller  than  a  \  it  is  crucial  to  take  this  temperature  dependent  a  for 
c-Si  into  account.  In  his  calculations,  Bell  [  2  ]  assumed  a  constant  /alue 
of  100/cm  for  absorption  coefficient  of  the  crystalline  substrate.  Such 
an  average  value  for  a  does  not  reflect  the  dynamics  of  the  processes  arising 
from  interband  and  free-carrier  absorption  processes.  However,  when  the 
amorphous  layer  is  thick  enough  (>  a  ^) ,  one  can  neglect  the  substrate 
effects  without  much  loss  of  accuracy.  The  absorption  coefficient  of  liquid 
Si  at  1.06  ym  is  assumed  to  be  1.4  x  10*V cm 

In  Fig.  3.1  we  have  plotted  the  temperature  at  the  surface  and  at 

2 

various  distances  inside  the  Si  at  different  times  for  a  fluence  of  2.5  J/cm  . 
At  10  ns,  the  temperature  profile  is  very  flat  indicating  that  the  tempera¬ 
ture  rise  is  being  governed  by  adiabatic  processes.  The  temperature  profile 
becomes  sharper  as  time  increases,  as  shown  by  the  profile  at  20  ns,  as 
thermal  diffusion  processes  become  more  important.  At  37.5  ns,  the  surface 
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temperature  rises  to  its  maximum  value  of  1780  K  and  the  gradient  in  the 
temperature  profile  is  maximum.  At  later  times  (e.g.,  80  ns)  the  profile 
becomes  very  flat  as  the  temperature  remains  constant  at  the  melting  point 
of  a-Si  to  a  depth  of  -1400  A. 

In  Fig.  3.2  we  show  the  time  evolution  of  the  melt  front  (x  )  for 

m 

three  different  laser  fluences.  We  note  that  as  the  laser  fluence  is 

increased,  the  initiation  of  the  melting  occurs  earlier,  as  expected. 

Also,  the  melt  front  propagates  deeper  into  the  material  with  increase  of 

2 

laser  fluence.  In  particular,  with  a  laser  fluence  of  3.0  J/cm  ,  the  melt 

front  propagates  to  a  maximum  depth  of  2800  A  Into  the  material,  which  is 

more  than  the  thickness  (1325  A)  of  the  Implanted  amorphous  region. 

From  Fig.  3.2  we  can  estimate  the  velocity  of  propagation  (v  )  of  the 

m 

melt  front  into  the  material  and  also  the  velocity  of  crystallization  (v^) 

as  the  melt  front  recedes  towards  the  surface.  We  find  that  for  values  of 

2 

fluence  of  2.0,  2.5,  and  3.0  J/cm  ,  the  value  of  v  is  535,  600,  and  950 

m 

cm/sec  respectively.  We  note  that  with  the  increase  of  the  laser  fluence, 

the  velocity  of  the  melt  front  also  increases.  The  velocity  of  crystalliz- 

2 

atlon  (v  )  from  the  laser  fluences  of  2.0,  2.5,  and  3.0  J/cm  are  found  to 
c 

be  90,  95,  and  100  cm/sec,  respectively. 

In  Fig.  3.3  we  plot  the  duration  of  the  surface  melting  and  the 
maximum  depth  of  the  melt  front  as  a  function  of  the  laser  fluence. 

From  time-dependent  reflectivity  measurements  we  found  that  the 

2 

threshold  for  melting  is  1.0  J/cm  .  We  have  adjusted  the  value  of  the 

absorption  coefficient  of  the  implanted  layer  such  that  the  calculated 

value  of  the  threshold  laser  fluence  for  melting  is  also  1.0  J/cm .  From 

Fig.  3.3  we  find  that  the  laser  fluence  necessary  to  melt  1325  A  is  2.3 
2 

j/cm  . 
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Figure  3.2.  Position  of  the  melt  front  vs.  time  for  different  laser 
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Figure  3.3. 


Calculated  values  of  the  duration  of  surface  melting 
and  the  maximum  depth  of  the  melt  front  vs.  laser 
fluence . 


Wafers  of  4  0-cm  phosphorus-doped  <100>  silicon  were  used  in  this 

investigation.  BF^  was  implanted  at  room  temperature  at  150  keV  to  a  dose 
15  2 

of  1  x  10  /cm  .  To  reduce  channeling  effects,  the  implantations  were 
performed  a  few  degrees  (  ~7°)  away  from  the  surface  normal.  This  implan¬ 
tation  produced  an  amorphous  layer  of  thickness  ~  1325  A,  as  determined  by 
HF  etching  experiments  [  31]. 

3.3  Time  Resolved  Reflectivity  Measurements 

The  details  of  the  experimental  set  up  for  doing  time-resolved  reflec¬ 
tivity  measurements  has  been  described  in  Chapter  2.3.  In  this  case  we  used 
a  5W  Ar  laser  operating  at  488  nm  and  a  5  mW  He-Ne  laser  operating  at  633  nm 
incident  on  the  sample  to  study  the  change  in  reflectance  during  the  annealing 
transient.  A  400  MHz  Tektronix  transient  digitizer  was  used  to  record  the 
transient  signals. 

A  typical  reflectivity  change  at  the  surface  probed  by  an  Ar  laser 

2 

(488  nm)  during  a  4.0  J/cm  Nd: glass  laser  pulse  is  shown  by  the  solid  line 
in  Fig.  3.4.  The  laser  pulse  is  Indicated  by  the  dashed  line.  We  see  from 
Fig.  3.4  that  the  reflectivity  increases  abruptly  from  its  value  for 
amorphous  silicon  at  room  temperature  ( ~ 40%)  to  a  value  of  66%  correspon¬ 
ding  to  the  value  of  liquid  silicon  at  488  nm  [  16  ].  The  value  of  reflec¬ 
tivity  remains  at  66%  during  the  molten  period  (180  ns) ,  and  then  decreases 
to  a  value  of  35%  corresponding  to  crystalline  Si  at  room  temperature. 

For  various  laser  fluences,  and  using  the  Ar  and  He-Ne  probe  lasers,  we 

measured  the  duration  of  the  flat  portion  (r  )  of  the  reflectivity  curve, 

m 

which  gives  the  time  that  the  surface  remains  molten.  Fig.  3.5  shows  the 


Nd:  Glass  Laser 
4.0  J/cm2 


Figure  3.4.  Reflectivity  profile  (solid  line)  of  the  Si  surface 

at  433  nm  as  it  passes  through  the  phase  transition  for 

2 

a  Nd:glass  laser  fluence  of  4.0  J/cm  .  The  dashed  line 
shows  the  laser  pulse  (27.5  ns  FWHM) . 
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plot  of  experimentally  determined  values  of  x  vs.  laser  fluence. 

There  have  been  speculations  by  Van  Vechten  et  al.  [7]  that  the 
phenomenon  of  pulsed  laser  annealing  is  a  nonthermal  process  resulting  from 
the  presence  of  a  dense  plasma  of  photoexcited  carriers.  As  discussed  by 
Nathan  et  al.  [9]  the  plasma  model  would  predict  a  probe  wave-length  de- 

4 

pendence  of  X  for  the  measured  melt-duration  T  .  However,  we  find  that  the 

measured  duration  of  surface  melting  x  is  independent  of  the  wave-length 

m 

of  the  probe  laser.  The  plasma  model  prediction  can  be  summarized  by  the 

relation  t  (488  nm)  ■  0.35  T  (633  nm) .  Since  our  experimental  results 
m  m 

conclusively  show  that  x^  is  independent  of  wavelength  as  probed  by  the  Ar 
or  He-Ne  lasers,  the  melting  model  properly  explains  the  dynamics  of  surface 
melting  and  crystalline  regrowth  during  the  pulsed  laser  annealing  studied 
here. 

The  melting  threshold  was  experimentally  determined  to  be  larger  than 
2 

1.0  J/cm  .  Samples  were  annealed  with  laser  fluences  of  1.0,  2.0,  2.5, 

2 

3.54,  and  4.0  J/cm  for  study  with  the  SIMS  and  the  Hall  effect  measurements. 

3. 4  Atomic  Profiling 

The  boron,  fluorine  atomic  profiles  presented  in  this  work  were  meas¬ 
ured  by  secondary  ion  mass  spectrometry  (SIMS)  I  41,  42]  using  a  CAMECA  IMS-37 
ion  microprobe  located  in  the  Materials  Research  Laboratory.  In  this  meas¬ 
urement  technique,  an  0^  primary  beam  was  rastered  over  an  area  of  ~250 
2 

urn  .  Standard  precautions  were  taken  to  reject  spurious  signals  from  the 

11  19 

crater  walls.  The  B  and  F  sputtered  ion  intensities  were  normalized 

30 

to  the  intensity  of  Si  sputtered  ions  from  the  matrix  to  compensate  for 
changes  in  the  secondary  ion  interaction  efficiency  as  different  specimens 


were  moved  into  the  analysis  position  of  the  ion  microscope.  Depth  scales 
were  determined  by  crater  depth  measurements  using  a  Dektak  stylus.  The 

concentration  scale  was  obtained  from  a  standard  prepared  by  Implanting 

+  15  2 

BF2  into  Si  to  a  dose  of  1  x  10  /cm  at  150  keV,  and  measuring  the 

unannealed  profile. 

Figure  3.6  shows  a  plot  of  the  boron  atomic  concentration  as  deter¬ 
mined  by  SIMS  measurements  for  the  as-implanted  sample  and  for  samples 

2 

annealed  with  laser  fluences  of  1.0,  2.0,  2.5,  3.5  and  4.0  j/cm  .  For  the 

2 

laser  fluence  of  1.0  j/cm  »  no  melting  was  observed  by  time-resolved 

2 

reflectivity  measurements.  The  boron  atomic  profile  for  1.0  J/cm  laser 

fluence  is  identical  to  that  of  the  as-implanted  sample.  Since  there  is 

2 

no  melting  at  1.0  J/cm  ,  there  is  no  redistribution  of  the  implanted 

profile,  due  to  the  small  diffusion  coefficient  of  boron  in  the  solid 

.  ,  . n-14  2  -lx 

phase  (-10  cm  s  ). 

2 

For  laser  fluences  of  2.0,  2.5,  3.5  and  4.0  J/cm  there  is  an  in¬ 
creasing  redistribution  of  the  boron  atomic  profile  with  Increase  of  the 
laser  fluence.  The  substitutional  redistribution  of  boron  induced  by  pulsed 
laser  annealing  as  shown  by  Figure  3.6  cannot  be  explained  by  thermal  dif¬ 
fusion  in  the  solid  because  the  time  scale  is  too  short.  Since  the  melt 
duration  is  on  the  order  of  a  few  hundred  nanoseconds,  this  redistribution 
of  the  boron  atomic  profile  can  only  be  attributed  to  diffusion  in  the  liquid 

phase  [4].  The  diffusion  coefficient  of  boron  in  the  liquid  phase  is 
4  2-1 

.  10  cm  s  .  As  the  laser  fluence  is  increased  beyond  the  melting  threshold 

2 

of  1.0  J/cm  ,  boron  diffuses  deeper  into  the  material  and  the  profile 
becomes  flatter  as  shown  in  Figure  3.6.  The  SIMS  results  shown  in  Fig.  3.6 
directly  support  the  thermal  melting  model. 

Figure  3.7  shows  a  plot  of  the  fluorine  atomic  concentration  as 


F19  Concentration  (cm-3) 
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determined  by  SIMS  measurements  for  the  as-implanted  sample  and  for  samples 

2 

annealed  with  laser  fluences  of  1.0,  2.5,  3.5  and  4.0  J/cm  .  The  profiles 

2 

for  the  as-implanted  and  the  1.0  J /cm~  samples  are  identical,  since  there 
is  no  melting.  The  most  significant  changes  in  the  fluorine  profiles  are 
the  appearance  of  pronounced  shoulders  in  the  profiles  at  -  2250  A  (for 

2.5  J/cm2),  -  2500  A  (for  3.5  J/cm2)  and  «  2700  A  (for  4.0  J/cm2). 

Apparent  changes  in  the  fluorine  concentration  near  the  peak  of  the  profile 
may  be  within  the  experimental  uncertainties  of  the  SIMS  technique.  The 
redistribution  of  fluorine  will  be  considered  below  along  with  the  electrical 
activation  of  the  boron  atoms. 

3. 5  Electrical  Profiling 

To  investigate  the  electrication  activation  after  pulsed  laser 
annealing  we  determine  the  sheet  resistivity,  sheet  carrier  concentration 
and  mobility  of  the  surface  layer.  Sheet  carrier  concentrations  are 
obtained  from  Hall  effect  measurements  using  a  van  der  Pauw  configuration. 
Differential  resistivity  and  Hall  effect  measurements  in  conjunction  with 
successive  layer  removal  are  used  to  obtain  the  electrical  carrier  distri¬ 
butions  as  a  function  of  the  laser  fluence.  We  observe  that  for  good 
quality  crystalline  regrowth  and  full  electrical  activation,  the  laser 
fluence  should  lie  between  the  melting  threshold  and  damage  threshold. 

3.5.1  Hall  System 

The  schematic  diagram  of  the  double  a.c.  Hall  system  [43]  is  shown 
in  Fig.  3.8.  A  sample  current  of  about  0.3  mA,  at  a  frequency  of  1  kHz 
and  a  200  Gauss  magnetic  field  at  250  Hz  were  used.  The  Hall  voltage 


Hall 

Coefficient 


Figure  3.8.  Block  diagram  of  the  double  ac  Hall  system. 


appears  at  the  heterodyne  frequency  of  750  Hz.  Buffer  stages  are  placed 
between  the  mixer  and  both  oscillator  stages  and  also  before  the  current 
Inputs  to  the  sample,  thereby  achieving  significant  noise  reduction.  Three 
stages  of  tuned  active  filters  are  used  to  generate  the  750  Hz  reference. 
The  final  two  stages  are  the  reference  channels  of  the  PAR  JB-5  and  HR-8 
lock-in  amplifiers. 


3.5.2  Sample  Preparation 

After  implantation  and  laser  annealing  the  samples  were  subjected 
to  anodic  oxidation.  An  anomalous  surface  layer,  sometimes  observed  in 
heavily  implanted  samples,  causes  a  rough  etched  surface.  This  surface 
layer  was  removed  by  anodic  oxidation  [44]  prior  to  chemical  etching.  In 
this  process  -  100  A  of  the  Si  surface  was  removed,  as  estimated  from 
measurement  of  the  anodic  oxide  thickness.  This  surface  treatment  was 
performed  on  all  samples  used  in  electrical  measurements  after  laser  anneal¬ 
ing. 

After  implantation,  laser  annealing,  and  anodic  oxidation  stripping, 
electrical  contacts  were  formed  by  evaporating  Au  through  a  shadow  mask 
and  then  sintering  at  200°C  in  flowing  ^  for  15  sec.  The  samples  were 
then  defined  photolithographically  in  a  van  der  Pauw  geometry.  Using  a 
diluted  planar  chemical  etch  [  45  ],  successive  Si  layers  were  removed.  To 
minimize  the  uncertainty  due  to  possible  variation  in  the  etch  rate, 
etched  step  heights  were  measured  at  various  points  in  the  profile,  using 
a  Dektak  stylus.  In  the  case  of  such  high  dose  B  implanted  Si,  the  ratio 
of  Hall  mobility  to  conductivity  mobility  was  taken  to  be  0.73,  appropriate 
for  high  concentration  p-regions  [46].  The  van  der  Pauw  geometry  and  the 
sample  holder  are  shown  in  Figure  3.9. 


Cross  sectional  view  of  the  sample  holder  used  for  double 
ac  van  der  Pauw  measurements.  The  insert  shows  the  sample 
configuration  and  spring  loaded  contacts  in  the  holder. 
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3.5.3  Data  Analysis 


The  sheet  resistivity  p  and  sheet  Hall  coefficient  R  are  given  by 

s  s 


V  +  V 
ABCD  BCDA 


./  VABCD  \ 
\VBCDA  I 


(3.3) 


B(rms)  I(rms) 


(3.4) 


Here  represents  the  voltage  between  contacts  A  and  B  with  current  I 

passed  between  C  and  D.  VBCDA  and  4VRnAr  are  similarly  defined,  with 


BDAC 


AVBDAC*  ttie  Hal1  volta8e  111  a  field,  B.  The  correction  factor 

f (VAn„,./V_„. )  is  Van  der  Pauw's  factor.  Let  the  sheet  resistivities 
ABCD  BCDA 


before  and 


after  stripping  a  layer  of  thickness  d.  be  donoted  by  p  and 

j  8j 


p  :  similarly,  R„  and  R_  refer  to  sheet  Hall  coefficients.  The 

sj-l  Sj  S3-l 

average  mobility  of  the  layer  is  given  by 


i  ,  1-1  "J-l  1 

P  P  P 

sj  8j-l  sj-l  '  sj 


(3.5) 


and  the  average  carrier  concentration  is  given  by 


P  P 

Jj _ i 


(3.6) 


Here,  the  Hall  coefficient  factor,  which  is  the  ratio  of  the  Hall  mobility 
to  conductivity  mobility,  is  assumed  to  be  0.73. 
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3. 5. A  Experimental  Results 

In  Figure  3.10  the  acceptor  concentration  and  mobility  profiles  are 

compared  with  the  boron  atomic  profile  for  the  sample  irradiated  with  a 

2 

laser  fluence  of  1.0  J/cm  .  We  observe  that  for  a  laser  fluence  of  1.0 
2 

J/cm  ,  there  is  nearly  full  electrical  activation  to  a  depth  of  'v  1300 

A,  the  boundary  of  the  original  amorphous-crystalline  Interface  (I) . 

Beyond  'v  1300  A,  there  is  no  electrical  activation.  We  also  note  from 

Figure  3.7  that  the  peak  fluorine  concentration  occurs  near  *v>  1200  A. 

2 

For  a  pulse  of  1.0  J/cm  ,  the  temperature  rises  to  a  value  below  the 

melting  point  of  amorphous  silicon,  producing  annealing  by  solid  phase 

epitaxial  regrowth.  The  termination  of  the  active  acceptor  profile  at  I 

is  consistent  with  results  of  furnace  annealing  [31]. 

Figure  3.11  -  3.14  shows  plots  of  the  acceptors  and  mobility  profiles 

2 

for  laser  fluences  of  2.0,  2.5,  3.5  and  4.0  J/cm  respectively.  All  these 
fluences  induce  melting  as  observed  by  time  resolved  reflectivity  measure¬ 
ment  s . 

Figure  3.12  shows  a  plot  of  the  electrical  carrier  profile  for  a 

2 

laser  fluence  of  2.5  J/cm  ,  which  produces  melting  and  redistribution  of 
the  boron  profile.  In  this  case  there  is  less  than  100%  electrical  acti¬ 
vation  to  a  depth  of  ^  2000  A.  The  peak  electrical  carrier  concentration 
is  a  factor  of  1.5  less  than  the  peak  atomic  profile.  We  note  from  Figure 
3.7  that  the  peak  fluorine  concentration  at  'v  1200  A  is  slightly  higher 
than  that  of  the  as-implanted  sample  and  there  is  a  second  smaller  fluorine 
peak  at  'v  2250  A,  beyond  which  we  find  no  electrical  activation.  Thus, 
there  appears  to  be  a  relation  between  electrical  activation  and  the 


fluorine  distribution. 


Electrical 

Carrier 

Profile 


B11  Atomic 
Profile 


V 


.Mobility 


<iOO>Si 
BF£,  150keV 
lxlOtycm2  R.T. 

Nd:  Glass  Laser  1.06^.m 
2.0  J/cm?  2.75ns 


Depth  (A°) 


UP*  1919 


Figure  3.11.  Net  acceptor  concentration  and  mobility  profiles  for 
the  sample  irradiated  with  2.0  J/cm^. 
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Figure  3.13  and  3.14  shows  similar  plots  of  acceptors  and  mobility 

2 

profiles  for  laser  fluence  of  3.5  and  4.0  J/cm  respectively.  For  this 
rather  high  laser  fluence  there  is  nearly  100%  electrical  activation  to 
a  depth  of  ^  2250  and  ^  2500  A  respectively.  Also  there  is  no  electrical 
activation  beyond  'v  2400  and  ^  2750  A  for  3.5  and  4.0  J/cm"  laser  flu- 
ences.  From  Figure  3.7  we  find  there  is  a  second  fluorine  peak  at  'v  2500 
A  (for  3.5  J/cm^)  and  ^  2700  A  (for  4.0  J/cm^). 

For  full  electrical  activation  with  pulsed  laser  annealing,  it  is 
absolutely  necessary  to  choose  the  laser  fluence  such  that  the  melt  front 
penetrates  deeper  than  the  defect  tail  beyond  the  amorphous-to-crystalline 
interface  [  47  ]. 

The  fluorine  concentration  after  the  laser  pulse  seems  to  be  influenced 
by  the  distribution  of  the  residual  defect  concentration.  The  above 
discussion  of  fluorine  migration  and  electrical  activation  profiles  clearly 
suggests  the  influence  of  residual  defects.  The  preferential  migration  of 
fluorine  to  a  defect -rich  region  is  similar  to  the  ion-damage  gettering 
discussed  by  Seidel  et  al.[48]  and  by  Tsai  e£  a_l.  [32].  For  all  the  laser 
fluences  the  damaged  region  responsible  for  fluorine  gettering  also  gives 
rise  to  the  electrically  inactive  tail  of  the  boron  distribution.  For 
the  2.5,  3.5  and  4.0  J/cm  fluences,  a  defect-rich  region  occurs  well  below 
the  original  amorphous-to-crystalline  interface  [  49  ]. 

3. 6  Electron  Channeling  Pattern  using  a  Scanning  Electron  Microscope 

To  investigate  the  material  properties  of  the  amorphized  layer 
before  and  after  the  laser  pulse,  we  employed  the  technique  of  electron 
channeling  using  a  scanning  electron  microscope.  For  the  as-implanted 
material  there  was  no  electron  channeling  pattern,  confirming  the  presence 
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of  an  amorphous  layer.  The  sample  that  was  irradiated  with  a  4.0  J/cm 
laser  fluence  clearly  showed  an  electron  channeling  pattern  due  to  the 
crystalline  regrowth  of  the  amorphous  layer.  The  regrown  layer  has  a 
<100>  orientation,  as  expected  for  recrystallization  on  the  < 100>  sub¬ 
strate.  A  typical  picture  of  the  channeling  pattern  is  shown  in  Fig. 

3.15. 

The  channeling  pattern  shown  in  Fig.  3.15  is  as  sharp  as  that  obtained 

from  an  unimplanted  sample  from  the  same  wafer.  This  supports  the  fact 

2 

that  a  fluence  of  4.0  J/cm  is  sufficient  for  the  melt-front  to  move  through 
the  original  amorphous-crystalline  interface,  as  predicted  from  the 


calculations. 
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4.  DYNAMICS  OF  PULSED  CO-  LASER 
ANNEALING  OF  SILICON  1 

Pulsed  lasers  have  been  successfully  applied  for  annealing  of 
radiation  damage  associated  with  ion-implantation  in  Si  during  device  fabrica¬ 
tion.  Experimentally,  work  has  been  mainly  confined  to  laser  wavelengths  * 
in  the  visible  and  near-infrared  using  a  frequency  doubled  YAG  laser  at 
X  *  0.  53  pm  [3  ]  and  a  Nd: glass  laser  at  X  =  1.06  pm  [36]  .  In  all  these 
lasers  the  spot-size  of  the  laser  beam  on  the  Si  sample  was  typically  a  few 
mm  in  diameter.  Hence,  to  anneal  a  large  area  it  was  necessary  to  use 
successive  overlapping  pulses,  leading  to  nonuniform  crystalline  quality  in 
the  overlap  regions.  Pulsed  CO2  lasers  have  the  distinct  advantage  of  having 
considerably  larger  spot-size,  and  are  therefore  very  promising  for  laser 
annealing  of  large  Si  wafers  commonly  used  in  industrial  production.  At 
10.6  pm,  the  absorption  coefficient  of  amorphous  silicon  is  much  smaller 
than  that  of  crystalline  silicon  having  large  donor  concentration.  Hence, 
for  annealing  with  a  CO^  laser  it  is  advantageous  to  implant  the  sample  with 
a  moderately  high  ion  dose  without  driving  the  sample  amorphous. 

We  present  a  detailed  calculation  using  the  method  of  finite  dif¬ 
ferences  to  study  the  dynamics  of  pulsed  CO2  laser  annealing  of  ion- implanted 
Si.  The  calculations  are  based  on  a  thermal  melting  model,  taking  into 
account  the  temperature  dependences  of  all  pertinent  material  parameters, 
including  the  absorption  coefficient.  We  investigate  the  "thermal  runaway" 
behaviour  and  calculate  the  threshold  power  density  to  cause  thermal  runaway 
and  surface  melting  of  Si  during  CO2  laser  annealing. 

We  present  results  from  a  quantum  mechanical  calcuation  of  the 
temperature  variation  of  the  free  carrier  absorption  coefficient  for  n-type 
Si  at  10.6  pm.  We  report  on  the  application  of  che  temperature  dependent 


absorption  coefficient  to  investigate  the  rise  in  surface  temperature 
leading  to  thermal  runaway  during  CO^  laser  annealing  of  Si.  We  have  cal¬ 
culated  the  power  density  required  to  initiate  the  thermal  runaway 
situation.  Our  results  are  in  good  agreement  with  those  observed  experi¬ 
mentally  [50,51].  In  calculating  the  absorption  coefficient  we  have  ad¬ 
justed  the  deformation  potentials  such  that  the  calculated  values  agree 
with  the  available  experimentally  obtained  values  [52,53]  at  300°  K.  It 
is  well  known  that  at  10.6  pm  the  absorption  mechanism  in  Si  is  primarily 
due  to  free  carriers,  with  a  small  contribution  due  to  multiphonon 
excitations.  The  variation  of  absorption  coefficient  with  temperature 
provides  insight  into  the  mechanism  by  which  laser  power  couples  into  the 
material.  We  find  that  the  absorption  coefficient  increases  rapidly  with 
temperature,  resulting  in  a  "thermal  runaway"  situation  for  large  power 
densities.  This  can  lead  to  melting  of  the  implanted  layer  as  observed 


experimentally,  and  subsequent  liquid-phase  epitaxial  regrowth  [51]  during 
pulsed  CO2  laser  annealing. 

We  calculate  the  contributions  from  the  acoustic  deformation 
potential  scattering  (a  ),  optical  deformation  potential  scattering  (a  ), 

2C  Op 

and  ionized  impurity  scattering  (0^) >  and  add  them  together  to  give  the 
total  free  carrier  absorption  coefficient  (a)  as  a  function  of  temperature 
from  300°  K  to  1685°  K  (the  melting  point  of  Si)  : 


a  *  a  +  a  +  a. 
ac  op  10 


(4.1) 


We  have  taken  into  consideration  the  contribution  from  both  the  electrons 


and  holes  for  completeness,  although  the  holes  contribute  little  for  n-type 
Si.  The  data  presented  are  for  the  case  of  thermal  equilibrium  with  the 


electron  temperature  equal  to  the  lattice  temperature  (T  *  T) . 
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The  free  carrier  absorption  coefficient  for  the  acoustic  deforma¬ 
tion  potential  scattering  is  given  by 


3/2  2  2  ,*  J-/ 2  r  -,1/ 2 

y  c  2  nee  (m  k_T)  _ 

o  ac  o  T  .  , 

- j— -  —  sinh 

/<  3n  h  c^o)  *■  e- 


"*BTe 


(4.2) 


where  y  c  ■  377  ohm,  n  is  the  free  carrier  concentration,  e  is  the  acoustic 
o  ac 

deformation  potential  constant,  <  is  the  dielectric  constant  of  the  material, 

* 

c ^  is  the  longitudinal  elastic  constant,  m  is  the  effective  mass  of 

electrons  (holes) ,  and  is  the  modified  Bessel  function  of  the  second 

order.  The  contribution  of  electrons  and  holes  toward  a  is  calculated 

ac 

by  using  the  parameters  given  in  reference  [54 ]  .  A  factor  of  2  is  taken 
into  account  in  equation  (4.2)  to  consider  both  the  phonon  absorption  and 
emission  processes. 

We  note  that  in  the  case  of  thermal  equilibrium  (T^  ■  T)  and  high 
temperature,  Eq.  4.2  reduces  to  that  obtained  by  using  the  classical  Drude 
model.  The  dependence  of  a  on  the  carrier  temperature  T  in  this  approxima- 

SO  6 

tion  is  given  by 


a  «  T 
ac  e 


(4.3) 


In  the  quantum  limit  hu>  »  2kgTe  the  absorption  becomes  indepen¬ 
dent  of  T  .  The  absorption  is  stronger  in  the  quantum  limit  than  it  is  in 

3 /2 

the  classical  limit,  but  the  dependence  on  X  is  less  strong  (*  \  ,  com¬ 


pared  with  A  in  the  classical  limit). 


The  total  concentration  of  electrons  is  obtained  by  assuming  that 

all  the  donors  are  Ionized  at  300°  K,  each  contributing  one  electron.  The 

number  of  holes  is  obtained  from  the  charge  neutrality  condition  and  the 
2 

relation  np  *  n^  .  The  total  free  carrier  concentration  is  given  by  the 

sum  of  electron  and  hole  concentrations. 

The  equations  for  the  optical  deformation  potential  scattering 

and  ionized  impurity  scattering  are  also  obtained  from  Ref.  55.  For  thermal 

2 

carriers  and  at  long  wavelengths,  a  *  X  .  This  is  the  same  as  for  acoustic 

op 

scattering.  For  a we  note  that  for  small  T^,  ctQ  is  nearly  independent 

1/2 

of  Tg,  whereas  for  large  Tg  a  is  proportional  to  T&  .  For  the  ionized 
impurity  scattering,  is  proportional  to  the  product  nN^,  where  Nj  is 
the  impurity  (donor)  concentration.  Since  for  an  uncompensated  semi¬ 
conductor  in  the  extrinsic  range  ■  n,  we  observe  an  increase  of  aio  with 
2 

n  rather  than  with  n  as  in  the  previous  cases  where  phonons  were  involved. 

Figure  4.1  shows  the  total  free  carrier  absorption  coefficient 
in  n-Si  as  a  function  of  carrier  concentration  for  various  wavelengths  at 
300°  K.  The  deformation  potentials  are  adjusted  so  that  the  calculated 
values  agree  with  the  experimentally  observed  values  [  52,53  ]  .  We  note  that 
for  a  particular  value  of  the  free  carrier  concentration  the  absorption  co¬ 
efficient  is  larger  for  longer  wavelengths.  Also,  for  a  particular  wave¬ 
length,  the  absorption  coefficient  increases  linearly  with  concentration. 

Figure  4.2  shows  the  variation  of  the  total  free  carrier  absorption 
coefficient  with  temperature  for  various  concentrations  of  ionized  donor 
impurities.  We  note  that  for  the  lower  values  of  the  donor  concentrations 
(Nj  -  10^/cm^,  10^/cra^)  and  T  <  700°  K,  the  total  number  of  free  carriers 

is  basically  given  by  the  number  of  electrons  contributed  by  the  ionized 


donors.  At  temperatures  T  >700  K,  however,  the  number  of  free  carriers 

contributed  by  the  intrinsic  carrier  concentration  becomes  comparable  to 

that  contributed  by  the  ionized  donors,  and  the  absorption  coefficient 

18  3 

rises  rather  sharply.  For  ■  10  /cm  ,  the  value  of  the  absorption 

coefficient  rises  from  a  value  of  70.0/cm  at  300°  K  to  a  value  of  4.0 
A  o 

x  10  /cm  at  1685  K.  This  rapid  increase  of  absorption  coefficient 
explains  the  thermal  runaway  observed  during  C02  laser  annealing  of  ion- 

implanted  silicon  [51]  .  We  also  note  that  at  1685°  K  the  values  of  a  for 

16  17  18  3  A 

=>  10  ,10  ,  and  10  /cm  approach  the  same  value  of  4.0  x  10  /cm,  , 

whereas  that  for  =  10^9 /cm3  approaches  a  value  of  5.0  x  loVcm-  So, 
for  efficient  coupling  of  the  laser  energy  to  the  material,  it  is  advisable 
to  start  with  Si  which  has  been  implanted  with  an  ionized  donor  impurity  of 
moderately  high  dose  (resulting  in  a  donor  concentration  of  1018  -  10^9/cm3) 
Alternatively,  the  sample  can  be  maintained  at  a  reasonably  high  temperature 
before  the  laser  energy  is  incident  so  that  enough  free  carriers  are  gen¬ 
erated  to  ensure  a  large  absorption  coefficient.  Thus,  Fig.  4.2  can  be 
used  in  choosing  appropriate  experimental  conditions  for  pulsed  or  cw  CO , 

laser  annealing  of  ion-implanted  Si.  The  results  of  Fig.  4.2  for  10.6  urn 

18  q 

can  be  fitted  by  an  empirical  power  law.  For  ■  10  /cm  ; 

a  -  1.016  x  10-2  T1,55,  300°  K  <  T  <  600°  K 

a  -  4.593  x  10’5  T2'38,  600°  K  <  T  <  800°  K 

a  -  5.559  x  10~19  T7,083,  900°  K  <  T  <  1685°  K 

A  similar  fit  can  be  obtained  for  the  other  donor  concentrations. 

These  approximations  can  be  used  in  choosing  the  laser  and  the  material 
parameters  and  substrate  temperature  for  successful  laser  annealing. 
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We  now  apply  our  calculations  of  the  temperature  dependent 
absorption  coefficient  to  Investigate  the  thermal  runaway  problem  during 
pulsed  CC>2  laser  annealing  of  Si. 

The  temperature  In  the  solid  is  obtained  from  the  solution  of 
the  one -dimensional  heat  flow  equation 


If  -  35  [*(T)f]+  E(1-R>- 


(4.4) 


where  ic  is  the  thermal  conductivity,  C  is  the  specific  heat  p  is  the 
density,  E  is  the  laser  power  density,  and  R  is  the  reflection  coefficient. 

For  simplicity,  we  assume  <  and  C  to  be  temperature  Independent. 


Equation  (4)  reduces  to 


3T  ^  32T  .  E  ( 1-R)  ae-CU£ 

at  -  3x2  cp 


(4.5) 


where  D  ■  ic/Cp. 


We  note  that  when  a  is  small,  i.e.,  the  absorption  length  is 
large  (a”^  >>  v'St) ,  we  can  neglect  the  heat  conduction  and  determine  the 
temperature  at  any  point  from  the  absorbed  energy.  This  is  no  longer  valid 
when  the  temperature  increases  to  the  point  that 

a-1(T  )  -  fit. 
c 

In  the  low  temeprature  approximation,  equation  (4.5)  simplifies 


3T  E ( 1-R)  -ox 

3t  "  Cp  06 


(4.6) 


We  now  apply  equation  (4.6)  to  find  the  temperature  at  the  sur¬ 
face  for  various  laser  power  densities  E. 

The  temperature  variation  of  the  absorption  coefficient  can  be 


expressed  by 
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T/T, 
a  ■  ae  1 


(4.7) 


■>19  . _ 3 


For  ■  10x  /cm  and  X  ■  10.6  pm,  we  find  a  *  317. 48 /cm  and  *  346°  K. 

The  solution  to  (6)  for  the  surface  temperature  is  then  given  by 


T(0,t)  *  Tq  +  T 


[a(Tn)E(l-R)t  ' 

1  pct^  ' 


(4.8) 


where  is  the  Initial  temperature.  For  Si  we  have  assumed  k  *  1.40  W/cmK, 

3 

C  *  0.74  J/cmK,  and  p  ■  2.33  g/cm  in  the  solid  phase.  The  value  of  R  at 
10.6  pm  is  0.30. 

To  investigate  the  thermal  runaway  behavior  we  have  plotted  in 

Figure  4.3.  The  surface  temperature  as  a  function  of  time  for  various  laser 

power  densities.  The  pulse  length  was  chosen  to  be  25  nanoseconds.  We 

observe  from  Figure  4.3  that  for  low  laser  power  densities  the  temperature 

ri.i*  is  not  fast  enough  to  initiate  thermal  runaway  during  the  25  ns  laser 

2 

pulse.  However,  for  power  densities  >75  MW/cm  we  observe  a  dramatic 

increase  of  the  surface  temperature,  leading  to  thermal  runaway.  This 

2 

minimum  power  density  E  *  75  MW/cm  is  in  good  agreement  with  observed 

experimental  values  [50,51].  We  note  that  as  the  absorption  coefficient 

becomes  large  with  Increased  temperature  we  cannot  apply  equation  (4.8).  In 

2 

Figure  4.3,  the  surface  temperatures  for  laser  power  densities  >  50  MW/cm 
have  been  calculated  up  to  a  value  T^  such  that  a”^(TC)  .  For  tempera¬ 

tures  higher  than  T^  equation  (4.8)  is  no  longer  valid  since  then  we  have 

the  condition  a  ^(T>T_)  <  /5t.  T_  is  given  by  the  expression  T.  *  T- in  — - — 1 

C  C  C  1  la/DtJ 

For  lower  values  of  N^,  the  absorptior  coefficient  is  smaller  and  the  minimum 
power  density  to  Induce  thermal  runaway  becomes  larger. 


Pulsed  C02  Laser 


± - 1 _ L 

5  10  15 

Time  (nanoseconds) 
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Surface  temperature  of  n-Si  with  *  10  /cm  during 
CO^  laser  pulses  (25  nsecs)  for  various  laser  power 


densities. 
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The  above  simplified  model  enables  us  to  calculate  only  the  sur¬ 
face  temperature  to  investigate  the  "thermal  runaway"  phenomenon.  We  now 
apply  the  melting  model  taking  into  account  the  temperature  variation  of  all 
the  material  parameters  to  Investigate  the  phenomenon  of  melting  and  regrowth 
with  a  pulsed  CC>2  laser. 

The  one-dimensional  heat  flow  is  described  by  the  equation 

C(T)  p  It =  h 1  k(t)  i£ ]  +  s(x»t)  (A-9) 

where  k  is  the  thermal  conductivity,  C  is  the  specific  heat  capacity,  and  p 
is  the  density  of  the  material.  The  source  term  S  is  governed  by  the 
mechanism  of  laser-solid  interaction  and  the  time  shape  of  the  laser  pulse. 
The  laser  pulse  is  assumed  to  be  Gaussian,  and  so  the  time-dependent  part 
of  the  source  term  is  given  by 

S(t)  -  4-  exp  (  -  (-^)2  ]  (4.10) 

ir”r  1 

where  T  is  the  parameter  which  determines  the  full  width  at  half  maximum 

(FWHM)  and  t'  is  the  time  when  the  pulse  reaches  its  maximum  value. 

The  spatial  part  of  the  source  term  is  given  by 

00 

S(x)  -  E  a(x)Cl-R)  exp  (-/  a(xf)dx’)  (4.11) 

o 

2 

where  E  is  the  laser  fluence  in  J/cm,R  in  the  reflectivity,  and  u  is  the 
temperature-dependent  absorption  coefficient.  During  melting  and  freezing 
the  latent  heat  L  has  been  accounted  for.  We  do  this  by  holding  the  tempera¬ 
ture  constant  until  enough  heat  is  absorbed  or  rejected  to  account  for  the 


latent  heat. 


To  solve  equation  (4.9)  numerically  by  the  method  of  finite  dif¬ 
ferences  we  divide  the  sample  into  equal  sections  of  length  Ax.  The  time 
step  is  At.  The  intervals  Ax  and  At  are  chosen  to  satisfy  the  von  Neumann 
stability  criterion  given  by 

At  <  (Ax)2 

specifically,  we  choose  Ax  *  250  8.  and  At  *  2  ps  for  our  calculations. 

The  temperature  dependence  of  the  free  carrier  absorption  co¬ 
efficient  at  10.6  pm  for  n-Si  has  been  calculated  [39]  by  taking  into  account 
acoustical  deformation  potential  scattering,  optical  deformation  potential 

scattering,  and  ionized  impurity  scattering.  For  an  ionized  donor  impurity 
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concentration  of  *  10  /cm  the  temperature  variation  of  the  absorption 
coefficient  is  given  by 

a(T)  -  «I/T 1 

2 

where  a  ■  3.17  x  10  /cm  and  T^  ■  346  K.  The  absorption  coefficient  in¬ 
creases  with  a  rise  in  the  concentration  of  carriers.  Also,  the  absorption 
coefficient  increases  nonlinearly  with  temperature,  leading  to  "thermal 
runaway"  and  thermal  melting. 
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A  typical  temperature  vs.  time  profile  for  n-Si  with  ■  10  /cm 
irradiated  with  a  pulsed  laser  of  25  ns  FWHM  and  laser  fluence  of  1.875 
J/cm  is  shown  in  Figure  4.4.  The  shape  of  the  laser  pulse  is  shown  in 
the  bottom  left-hand  corner.  The  notation  T^  denotes  the  melting  point  of 
silicon,  which  is  1685  K.  We  note  that  for  temperatures  T<700  K,  when  a  is 
small  and  the  absorption  length  is  large,  the  temperature  rises  rather  slowly 
and  is  mainly  determined  by  the  absorbed  energy  with  negligible  thermal 
conduction.  However,  as  ct  increases  rapidly  with  temperature,  the  slope  of 
the  temperature  rise  becomes  steeper,  heat  is  more  absorbed  in  the  near  surface 


Temperature 


2 

a  laser  fluence  of  1.875  J/cm  is  sufficient  to  initiate  "thermal  runaway", 
eventually  leading  to  surface  melting.  We  also  note  that  the  peak  of  the 
laser  pulse  occurs  at  25  ns,  whereas  the  maximum  surface  temperature  (1930  K) 
occurs  at  37.5  ns.  From  Figure  4.4  we  find  that  the  front  surface  remains 
molten  for  a  period  of  time  =  60  ns.  This  time  can  be  verified  experi¬ 
mentally  using  similar  time-resolved  reflectivity  measurements  as  discussed 
for  a  Nd:glass  laser  [36]  .  In  Figure  4.4  we  have  also  plotted  the  temperature 
profile  at  0.5  pm  at  1.0  pm.  We  note  that  deep  into  the  material  the  maximum 

temperature  is  smaller  and  the  temperature  peak  occurs  at  a  later  time,  as 

2 

expected.  For  1.875  J/cm  ,  we  find  from  our  data  that  the  melt  front  needs 
a  maximum  distance  of  0.15  pm. 

In  Figure  4.5  we  have  plotted  the  temperature  profile  in  the 

sample  at  different  times  for  a  laser  fluence  of  2.5  J/cm  .  The  maximum 

surface  temperature  of  2430  K  occurs  at  35  ns.  We  note  that  in  the  beginning 

of  the  pulse  the  slope  of  the  curve  is  rather  flat,  indicating  that  for  small 

a(large  absorption  length),  the  sample  is  heated  rather  uniformly.  However, 

as  the  temperature  increases,  the  absorption  length  decreases  and  the  surface 

is  heated  rapidly.  The  strongest  variation  of  temperature  with  distance 

occurs  at  30  ns,  when  the  surface  temperature  reaches  its  maximum  value. 

With  further  increase  in  time,  the  curve  becomes  flatter  again  as  heat  is 

transmitted  away  by  thermal  conduction. 

Figure  4.6  shows  a  plot  of  the  maximum  surface  temperature  and  the 

maximum  depth  of  melting  (x^)  for  various  laser  fluences.  We  note  that  a 

2 

laser  fluence  of  1.25  J/cm  causes  "thermal  runaway"  but  the  surface  tempera¬ 
ture  rises  to  a  maximum  value  of  only  920  K.  To  have  surface  melting,  a 
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minimum  laser  fluence  E  .  of  ~  1.52  j/cm  is  necessary.  As  the  laser 

min 

fluence  is  increases  beyond  E^^,  the  maximum  surface  temperature  increases 

and  the  depth  of  melting  also  increases.  Typically,  to  melt  a  layer  of  0.4 

ym  requires  a  laser  fluence  of  2.45  j/cnr,  as  can  be  seen  from  Figure  4.6. 

Figure  4.7  shows  the  variation  of  the  duration  of  surface  melting 

(t  )  as  the  laser  fluence  is  increased  from  the  threshold  value  for  melting 
m 

2  2 

(1.52  J/cm  ).  We  find  that  for  fluence  1.52  <  E  <  1.65  J/cm  the  value 

2 

of  increases  rather  sharply.  Beyond  1.65  J/cm  the  slope  of  the  curve 

2 

decreases.  However,  for  laser  fluences  >  1.85  J/cm  ,  the  slope  of  the  curve 

again  increases.  This  variation  of  t  is  comparable  to  the  variation  of 

m 

2 

the  maximum  surface  temperature  for  1.52  <E<  2.5  J/cm  as  shown  in  Fig. 

4.6.  Once  the  laser  fluence  increased  beyond  the  melting  threshold,  we 

would  expect  a  marked  variation  of  x^  with  E  as  shown  for  1.52  <  E  <  1. 65 
2 

J/cm  .  The  change  in  x^  for  the  three  different  regions  of  the  curve  is 

primarily  governed  by  the  magnitude  of  the  laser  fluence  relative  to  the 

threshold  value  (E  for  melting. 

Our  calculations  [56]  give  a  physical  picture  of  the  dynamics  of 

pulsed  CO^  laser  annealing.  The  model  provides  values  of  laser  fluences  or 

power  densities  required  to  melt  a  particular  thickness  for  a  laser  pulse 

19  3 

of  25  ns  (FWHM).  The  graphs  presented  here  are  for  ■  10  /cm  ;  one  can 
calculate  similarly  the  melting  threshold  E^^  for  other  concentrations  of 
ionized  donor  impurities.  We  note  that  since  a  decreases,  the  value  of  E 

min 

consequently  increases. 
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APPENDIX  A 


COMPUTER  PROGRAM  FOR  CALCULATION  OF  THE  DYNAMICS  OF  SURFACE 
MELTING  DURING  PULSED  LASER  ANNEALING  OF  AN  INFINITE  ION- 
IMPLANTED  AMORPHIZED  SEMICONDUCTOR 

This  program  is  titled  MELT. FOR.  This  program  is  very  versatile. 

It  can  be  used  for  calculation  of  the  temperature  rise  and  the  motion  of 

the  melt  front  during  pulsed  laser  annealing  of  a  semiconductor. 

For  a  particular  semiconductor,  e.g..  Si,  Ge,  GaAs,  one  has  to 

supply  data  for  the  material  parameters.  The  coupling  of  the  laser  and 

the  semiconductor  is  primarily  governed  by  the  absorption  coefficient. 

The  laser  pulse  width  is  characterized  by  GAMMA  and  the  f luence  is  given 

2 

by  energy  density  in  J/cm  . 

This  program  can  be  appropriately  modified  to  take  into  account  the 
substrate  effects.  This  situation  arises  when  one  considers  an  Implanted 
amorphized  layer  of  a  few  thousand  angstroms  thickness  on  a  crystalline 


substrate. 
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't 

'j 

0  CALCULATION  0?  THE  DYNAMICS  OP  SURFACE  MELTING 

0  DURING  LASER  ANNEALING  0?  AN  INFINITE  A-LAYER 

DIMENSION  T3NEW  (0:200), TSOLD (0 : 200 ) , X( 0 : 200 ) , CS ( 0 : 200 ) , 

1  AKS ( 300 : 2000 ) , ACS ( 300 : 2000 ) , 3KS ( 300 : 2000 ) , 

2  PI (300:2000) , P2 ( 300 : 2000 ) ,3CS( 300: 2000) , 

3  TLNEW( 0:200) ,TL0LD(0 : 200) ,TP(0:200) 

REAL  KS(0: 200) ,KL,KL1 , L, KLEFT ,KRIGHT 
n 

C  INPUT  DATA 

ft 

WRITS  (5,300) 

500  FORMAT  ('030LID  DENSITY  IN  GM/CC  ?") 

READ  (5,100)  RHOS 

ioo  format(ei i .3) 

n 

V 

WRITE  (5,505) 

505  FORMAT  ('OLIQUID  DENSITY  INGM/CC  ?') 

READ  (5,105)  RHOL 
105  format  (F1 1.3) 

-1 

'WRITE  (5,510) 

510  FORMAT  ( '0A3S0RPTI0N  SOLID  ?') 

READ  (5,* 10)  ALPHAS 
HO  FORMAT  (J10.2) 

n 

j 

WRITE  (5,515) 

515  FORMAT  ('0AB30RPTI0N  LIQUID  ?') 

READ  (5,115)  ALPHAL 
115  FORMAT  (310.2) 

r* 

■j 

WRITE  (5,520) 

520  FORMAT  ('OREFLECTION  SOLID  ?') 

READ  (5,120)  REFS 
120  FORMAT  (710.2) 

'"'l 

v> 

WRIT3(5,525) 

525  FORMAT ('OREFLECTION  LIQUID  ?') 

READ  (5,125)RSFL 
125  FORMAT (F1 1.3) 

C 

WRITE (5, 530) 

530  FORMAT  ('OENERGY  DENSITY  ?') 

READ  (5, 130) ENERGY 
130  FORMAT ( FI  0 . 2 ) 

r\ 


540 


WRITS  (5,540) 

FORMAT  ('OGAMMA  ?') 

READ  (5,140)  GAMMA 
HO  FORMAT  (BIO. 2) 

WRITE  (5,550) 

550  FORMAT  ('ODX  ?') 

READ  (5,150)  DX 

150  FORMAT  (S1D.2) 

"1 

WRITE  (5,550) 

560  FORMAT  ('ODT  ?') 

READ  (5,150)  DT 
160  FORMAT  (E10.2) 

C 

WRITE  (3,500)  RHOS.RHOL, ALPHAS, ALPHAL, REFS, REFL, 

2  ENERGY , GAMMA , DX , DT 

600  FORMAT ( 1  HI ,'RH0S=',F11 .3/1H0, 'RH0L-' ,F1 1 .3/ 

2  1H0, 'ALPHAS-' ,310. 2/lHO, 'ALPHAL-'  ,310.2/ 

3  '0REF3- ' , FI  0 . 2 / 'OREFL- ' , F 1 0 . 2 / 1  HO , 'ENERGY- ' , 

4  ? 10.2/1  HO, 'GAMMA*' ,E1 0.2/1  HO, 'DX-' ,E10.2/ 

5  ISO, 'DT*', BIO. 2) 

n 

C  CALCULATE  THERMAL  CONDUCTIVITIES  AND  SPECIFIC  HEATS 

q 

KL-0.322 

CL-0.915 

L-4220 

r* 

DO  50  1-300,2000 
IF( I.GT.600)  GO  TO  1 
PI  (I)-1336.35*I**(-1 .21 ) 

P2  ( I )-0. 19*I**0.24 
GO  TO  50 

1  I?(I.GT.900)  GO  TO  2 

PI  ( I)* 3770. 63*1** (-1  .36) 

P2(l)-0.19*I**0.24 
GO  TO  50 

2  IF ( I . GT . 1 200 )  GO  TO  3 
PI  ( l)- 1 031 .0*I**(-1 .15) 

?2(I)-0.41*I**0.12 

JO  TO  50 

5  ?i  (I't-’S.  n*:**(-o.53) 

?2(l)-0.  41  *I**0. 1  2 
50  CONTINUE 


00  51  1-3  00,1999 


ca  ca  o  -*•  ca  ca  ca  o  va  ca  ca  -»  ca  ca  ca  cn 


AKS(l)-P1 (1+1 )-?1 (I) 
ACS(l)*P2(l+1 )-P2(l) 
BKS(l )*P1 (l)-AKS(l)*I 
3CS(l)=P2(l)-ACS(l)*I 
CONTINUE 


INITIALIZATION 

DO  10  1=0,200 
TSOLD(I)=300. 

X(l)-FL0AT(l)*DX 
KS(l)-1 .395 
CS(l)-0.7469 
CONTINUE 

NMELT-0 

NFRONT-O 

M-1 

X30LD-Q.0 
XBNEVf =0 . 0 

CALCULATE  CONSTANTS  USED  IN  THE  CALCULATION 


DX1 -DX*2 . 

DTI -DT/RH0S 

PS1=0.554*(l .-EXP(-ALPHAS*DX))/DX*3NERGY*(l . -REFS ) / GAMMA 
DT  2  -DT  /  ( RH0  L*  CL ) 

DT3=L/DT 

PL1 -0.564*0  . -EXP( -ALFHAL*DX) )  /DX*ENERGY*(  1  .-REFL)  OAMMA 
KL1-KL/DX 

START  OF  THE  TIME  ITERATION 

DO  20  MM-50,500,50 
1  TIME-FLOAT (M)*DT 

TEXP-EXP( -( TIME /GAMMA- 1  .5S5)**2) 

IF(NMELT.EQ.O)  GO  TO  70 

CALCULATE  NEW  XB 

3 L0P33-KS ( N30LD) * ( TSOLD ( NBOLD ) -1  335 .)  -'(FLOAT (NB0LD)*DX-T30LD'> 
IF(X30LD.EQ.FL0AT(NB0LD-1 )#DX)  GO  TO  15 

3L0PEL-KL*(1635.-TL0LD(NB0LD-1 ))/(X30LD-FL0AT(N30LD-i  )*DX) 

GO  TO  13 

15  3L0PEL-KL*  ( TLOLD  ( NBOLD  - 1  )-?L0LD(NB0Lu-2)VDX 

1 3  X3NEW-X30LD+( 3L0PE3-3L0PSL) /DT3 

NBOLD- I?IX(X30LD/DX)+1 


N3NStf=IFIX(X3NEN/DX)+1 
IF(NFR0NT.EQ.1 )  50  TO  12 

CHECK  THE  SION  OP  FRONT  VELOCITY 

IP ( X3NEW . OE . X30LD)  00  TO  12 
VRITS(3,645)  TIME , X30LD 

F0RMAT(1H1 MELT  FRONT  REACHES  MAXIMUM  AT  IPS 11. 3, 
2  '  SEC'/IHO,'  X3MAX-  '.B11.3) 

NFR0NT*1 
00  TO  12 
NBOLDO 
N3NEWO 

TEMPERATURE  IN  THE  SOLID 


DO  30  N=NB0LD+1 , 1 39 
33=’SXP(-ALPHA3*(X(N)-.5*DX))*TSX? 

IF(N.EQ.I)  B3-TSXP 
P3-P31»33 

I?( NMELT.EQ. 1  )  PS-0.0 
XLEFT=*0  . 5*(XS(N-1  )+KS(N)) 

KRIGHT-0.5*(K3(N+1 )+K3(N)) 

A=-KLEFT*(TS0LD(N)-T30LD(N-1 )) 

2  +KRI0HT*(TS0LD(N-H  )-TSOLD(N)) 

TSNEW(N)=T30LD(N)  +DT1  ''C3(N)*(PS+A/(DX*DX) ) 

CONTINUE 
TSNEW ( 200 )  =  300 . 

IF(NMELT.EQ.O)  00  TO  33 

I?(N3M3W . OE  .TBOLD-H )  00  TO  5V 

DELT3M*(l  533. -T3N3N(  BOLD*1  )) /'?L0  VT(:l30LD+1  )-t3HEN/DT) 
DO  53  N=N3 OLD , N3NE.N ,  -> 

T3U3N  ( N ) = I3NEN  ( H  1  VDELTSM 
CONTINUE 

IF v \ NFRONT .EO • 0) . OR. (N3NE.N . OE .3 ) ^  00  TO  63 
SURFACE  MELTING  ENDS 
NR ITS ( 5 ,640)  TIME 

FORMAT ( I  HI SURFACE  MSLTINO  ENOS  AC  ',’?E'1.3,'  SEC'O 

X3NSW“0 . 0 

X30LD-0.0 

N3NEWO 

NMELT*0 
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N?R0NT=0 

DO  51  N=0 , N3NEW-1 
TSNS>f(M)=t635. 

00  TO  200 

I? ( T3NSW ( 1 ) . LT . 1 000 . )  00  TO  30 

TSNEW ( 0 ) -T3NEW ( 1 )-5 . 331 0»SXP( -4 • 3E4/T3NEW( 1 ) )*DX/KS( 1 ) 
GO  TO  200 
TSNEW(0)=T3NEW( 1  ) 

GO  TO  200 

TEMPERATURE  IN  THE  LIQUID 

DO  50  N=1 ,N30LD-2 

B4=EXP( -ALPH AL*( X ( N ) - . 5*DX ) ) *TEX? 

IF(N.EQ.l)  B4=TSX? 

PL=PL1 *B4 

A1 =KL1*(TL0LD(N+1 )-2.*TL0LD(N)+TL0LD(N-1  )) 
TLNEL(N)=TL0LD(N)+DT2*(?L+A1  /DX) 

CONTINUE 

I?(TLNE>T(1  ).LT.2000.)  GO  TO  31 

TLNEW ( 0 ) * TLNEW ( 1  )-5  •  3E1 0*SXP(  -4 .3E4/TLNEW(  1  )  )  /XL1 
GO  TO  32 

TLN3W(0)*TLNEW(l  ) 

IF(NBNEW. LS. (N30LD-1 ))  GO  TO  54 

DELTL=* ( TLNEW ( N30LD -2 )  -1  635  • )  /  (X3NEW/DX -FLOAT  (N30LD-2  ) ) 
DO  52  N=NB0LD-2 ,NBNEW-2 
TLNEW (N+1  )  =>  TLNEW  ( N )  -DELTL 
CONTINUE 

CALCULATE  NEW  THERMAL  CONDUCTIVITY  AND  SPECIFIC  HEAT 

DO  52  N=N3NSW,200 
INDEX“I?IX(T3NEW(N) ) 

IF ( INDEX- 300)  700,701,701 
KS ( N ) » -5 . O^E-3  *TSNEW ( N ) >2 . 347 
C3(N)«4. 52273-4 *TSNEW(N)+0. 51 3099 
GO  TO  52 

IF (INDEX- 2 000)  702,703,703 
•<3(N)--4. 3S-5*TSNEW(  N )  +0 .2469 
C3(N)-9.54533-5*T3NEW(N)0.35  91 
GO  TO  52 

X3 ( N ) ■ AK3 ( INDEX ) *TSNEW ( N ) + 3K3 ( INDEX ) 
CS(N)»AC3(INDEX)*T3NEW(N)+3C3( INDEX) 

CONTINUE 

IF(NM3LT.EQ.1 )  GO  TO  207 


CHECK  INITIAL  MELTING 


L 


1.25  Mil  1.4 


MICROCOPY  RESOLUTION  TEST  CHART 
’  NATIONAL  BUREAU  OF  STANDARD5-1963-A 


O  O  O  CJ  o  u  o 


IF(TSNBW(0)  .LT.  1 635  • )  GO  TO  205 
DO  71  N-0,200 

I? ( TSNEW ( N )  -1  635 • )  72,71,71 

71  CONTINUE 

72  XBNEW  »DX*FL0 AT ( N- 1  ) 

2  +DX*( TSNEW (N-1 )-1 635 . )/(TSNEW(N-1  )-TSNEW(N)) 

I?(X3NEW.LT.(2.*DX))  GO  TO  205 
WRITS(3,535)  TIME,X3NEW 

555  ?0RMAT( 1  HI  , '  3U3PAC3  MELTED  AT  '  ,  1  ?31 1  . 3 , '  330' 

2  /HO,'  XBNE Bf«  ',311.3) 

NMELT-1 

N3MSW-N 

DO  73  N“0,N3NBW-1 


73 

TLNBW (  M ) *TSNEW ( N ) 

REPLACE  OLD  VALUES 

BY  NEW  VALUES 

207 

XBOLD*XBNEW 

N30LD=*NBNEW 

DO  206  N=0,NBNEW-1 

205 

TL0LD( N ) =TLNEW ( N ) 

205 

DO  203  N-NBNEW.200 

203 

T50LD(N)*TSNEW(N) 

210 

M»M+1 

IF(M-MM)  11,40,40 

PRINT  OUT  TIME,  XB 

AND  TEMPERATURE 

40 

DO  41  N*NBNEW , 200 
TP(N)»TSOLD(N) 

41 

CONTINUE 

IP(NBNEW.SQ.O)  GO  ' 
DO  44  N*0 ,NBNEW-1 

TO  43 

44 

TP(N)-TLOLD(N) 

43  WRITE(5,510)  TIME 

610  F0RMAT(lH1 , 'TIME-' ,1  PE1 1 .3) 

WRITE(5 ,61 5)  XBNEW 
515  F0RMAT(lH0,'X3>',1PE11 .3) 

WRITE (5, 6 20) 

620  F0RMAT( 1 H0,4(5X , 'X' ,3X , 'TEMP' ,4X) ) 

DO  42  N-0,49 
N1-N+50 
N2-N1 +50 
N3-N2+50 
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